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PROGRAM  SUMMARY 


The  objective  of  the  program  is  to  develop  the  design  criteria  and 
analytical  methods  necessary  to  ensure  the  damage  tolerance  of  aircraft 
attachment  lugs.  As  planned,  the  program  proceeds  logically  from  an  exten¬ 
sive  cracking  data  survey  and  nondestructive  inspection  (NDI)  assessment, 
through  method  development  and  evaluation,  to  the  preparation  of  damage 
tolerance  design  criteria  for  aircraft  attachment  lugs. 

The  program  consists  of  three  phases  involving  seven  tasks.  Phase  I 
consists  of  Tasks  I,  II  and  III;  Phase  II  consists  of  Tasks  IV,  V  and  VI; 
and  Phase  III  consists  of  Task  VII.  A  roadmap  shown  in  Figure  l-l  summar¬ 
izes  the  major  activities  by  task,  decision  points  and  their  interrelation¬ 
ships. 

Task  I  involves  a  survey  of  structural  cracking  data  such  as  the  ini¬ 
tial  flaw  size,  shape  and  location  which  occur  in  aircraft  attachment  lugs. 
Sources  for  these  data  Include  open  literature,  available  Lockheed  data, 
and  visits  to  the  five  Air  Force  Air  Logistics  Centers  (ALCs).  The  types 
of  aircraft  structure  used  to  obtain  these  data  include  service  aircraft, 
full-scale  test  articles,  component  test  articles,  and  coupon  specimens. 

Task  II  assesses  the  current  NDI  capability  to  find  these  flaws  or 
cracks.  This  assessment  is  to  be  based  upon  information  obtained  from  the 
open  literature,  available  Lockheed  NDI  data  and  experience,  and  Air  Force 
ALC  data.  The  NDI  techniques  capable  of  finding  flaws  in  attachment  lugs 
and  the  flaw  sizes  these  techniques  are  capable  of  finding  are  identified. 
Where  possible,  the  probability  of  detecting  a  flaw  of  a  particular  size 
for  the  NDI  technique  involved  is  specified  as  well  as  the  confidence  level 
assigned  to  that  probability.  The  results  obtained  from  Tasks  1  and  11  will 
be  used  in  the  formulation  of  the  initial  flaw  assumptions  to  be  developed 
in  Task  VII  as  part  of  the  damage  tolerant  design  criteria  for  attachment 
lugs. 

Task  Itl  involves  three  different  levels  of  complexity  and  degrees 
of  sophistication  for  determining  stress  intensity  factors  for  single 
corner  cracks  and  single  through-the-thi ckness  cracks  in  aircraft 
attachment  lugs,  aud  the  development  ot  crack  growth  analysis  capable  of 
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predicting  the  growth  behavior  of  these  cracks  and  residual  strength  of 
these  lugs.  These  stress  intensity  factors  and  crack  growth  analyses  are 
used  in  Task  IV  to  predict  the  residual  strength  and  the  crack  growth  behav¬ 
ior  for  a  number  of  different  geometries  and  test  conditions  defined  in  the 
experimental  program.  These  predictions  are  made  prior  to  testing.  Two 
groups  of  attachment  lug  geometries  are  tested  and  experimental  test  data 
are  generated  in  Task  V.  The  analytical  methods  developed  in  Task  111  are 
evaluated  by  correlating  the  analytical  predictions  made  in  Task  IV  with 
the  Group  I  experimental  test  data  (based  upon  accuracy  and  cost)  for  use  in 
prediction  of  Group  II  tests.  Further  evaluation  of  the  selected  method  is 
made  by  correlating  the  analytical  predictions  for  the  Group  II  tests 
(Task  IV)  with  the  experimental  test  results  (Task  V).  These  correlations 
indicate  what  improvements  are  necessary  for  the  selected  analytical  method. 
The  results  are  presented  in  parametric  format  useful  to  designers  and 
analysts.  Damage  tolerant  design  criteria  for  aircraft  attachment  lugs  are 
developed  in  Task  VII.  These  criteria  are  similar  in  nature  to  those  of 
Military  Specification  MIL-A-83444,  and  require  crack  growth  analyses  by  the 
types  of  methods  developed  and  verified  in  Tasks  III  through  VI.  The  cri¬ 
teria  include  initial  flaw  assumptions  (e.g.,  initial  flaw  type,  shape, 
size,  etc.)  based  upon  the  cracking  data  survey  of  Task  I,  NDI  assessment 
of  Task  II,  and  crack  initiation  tests  of  Task  V. 

As  Figure  1-1  shows,  the  following  sequence  of  final  report  volumes  is 
generated  under  this  project: 

Volume  I.  Cracking  Data  Survey  and  NDI  Assessment  for  Attachment 
Lugs 

Volume  II.  Crack  Growth  Analysis  Methods  for  Attachment  Lugs 

Volume  III.  Experimental  Evaluation  of  Crack  Growth  Analysis  Methods 
for  Attachment  Lugs 

Volume  IV.  Tabulated  Test  Data  for  Attachment  Lugs 

Volume  V.  Executive  Summary  and  Damage  Tolerance  Criteria  Recommen¬ 
dations  for  Attachment  Lugs 

Volume  Vi.  User's  Manual  for  "LUGRO"  Computer  Program  to  Predict 
Crack  Growth  in  Attachment  Lugs 

This  is  Volume  IT  on  Crack  Growth  Analysis  Methods  for  Attachment  Lugs, 
which  is  the  result  of  Task  III  efforts.  It  contains  the  developed  analyt¬ 
ical  methodologies  for  determining  the  stress  intensity  factors  for  single 
corner  cracks  and  single  through-the-thickness  cracks  in  aircraft  attachment 
lugs,  and  for  predicting  the  growth  behavior  of  these  cracks. 
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SECTION  II 
INTRODUCTION 


In  aircraft  structures,  lug-type  joints  are  frequently  used  to  connect 
maj-'r  structural  components  or  in  linkage  structure.  The  lug  joint  is  norm¬ 
ally  connected  by  a  single  bolt  or  pin,  creating  a  simple  joint  that  is  easy 
to  assemble  and  disassemble.  Since  clamping  of  the  joint  is  not  normally 
allowed,  the  lug  can  act  as  a  pivot.  But  the  elastic  gross  section  stress 
concentration  for  normal  lugs  is  very  high,  resulting  in  a  relatively  short 
crack  initiation  and  crack  growth  life.  To  improve  the  crack  initiation 
life,  the  stress  concentration  factor  can  be  effectively  reduced  by  cold 
working  the  hole  or  by  installing  an  interference-fit  bushing  prior  to  pin 
fitting.  To  minimize  the  wear  in  either  the  lug  or  the  pin,  most  aircraft 
lugs  have  an  oil  fitting  and/or  lubrication  provisions. 

During  the  past  decade,  the  influence  of  fracture  mechanics  on  the 
design,  manufacture,  and  maintenance  of  aircraft  has  steadily  increased. 
Also,  nondestructive  inspection  techniques  have  been  improved  significantly. 
However,  some  cracks  still  cannot  be  detected  during  routine  maintenance 
inspection.  Under  service  loading,  such  cracks  will  grow  and  fracture  can 
occur  if  the  crack  length  reaches  a  critical  dimension  before  it  can  be 
detected  and  the  part  repaired  or  replaced.  To  assure  aircraft  safety,  the 
U.  S.  Air  Force  has  imposed  damage-tolerance  requirements  (MIL-A-83444 /[l ] 
which  include  the  prediction  of  fatigue  crack,  growth  life  and  residual 
strength  of  the  structure  by  assuming  that  small  initial  flaws  exist  at 
critical  locations  of  new  structure  due  to  various  material  and  manufactur¬ 
ing  and  process  operations.  Assumptions  regarding  the  initial  size,  shape, 
location,  multiplicity,  etc.  for  these  flaws  are  specified  in  MIL-A-83444. 
However,  these  assumptions  were  established  primarily  for  skin-stringer 
structure  and  may  not  be  applicable  for  attachment  lugs. 

Attachment  lugs  are  some  of  the  most  fracture  critical  components  in 
aircraft  structure,  and  the  consequences  of  a  structural  lug  failure  can  be 
very  severe.  Therefore,  it  is  necessary  to  develop  damage  tolerance  design 
requirements,  similar  to  MIL-A-83444,  for  attachment  lugs  to  ensure  the 
safety  of  aircraft.  The  development  of  these  damage  tolerance  requirements 
will  be  based  upon  actual  cracking  data  for  attachment  lugs  and  current 
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nondestructive  inspection  capability.  Once  the  damage  tolerance  design 
requirements  for  aircraft  attachment  lugs  are  established,  the  analytical 
methods  necessary  to  satisfy  the  crack  growth  and  residual  strength  require¬ 
ments  are  needed.  In  particular,  stress  intensity  factors  for  cracks  in 
attachment  lugs  are  needed.  Such  stress  intensity  factors  will  depend  upon 
the  complexities  of  structural  configuration,  crack  geometry,  applied  loads, 
and  the  fit  between  the  pin  and  the  lug. 

There  are  a  number  of  different  methods  for  determining  the  stress 
intensity  factors,  K,  for  cracks  in  aircraft  attachment  lugs.  Schijve  and 
Hoeymakers  [2]and  Wanhill [3]  derived  empirical  K-solutions  from  the  growth 
rate  data  for  through  cracks  under  constant  amplitude  loading  using  a  back¬ 
tracking  method  such  as  that  proposed  by  James  and  Anderson[4],  Analytic¬ 
ally,  Liu  and  Kan  [  5]and  Kirkby  and  Rooke  [6]  used  the  simple  compounded  solu¬ 
tion  method  which  involves  superimposing  known  solutions*  such  as  in  Refer¬ 
ence^]  to  estimate  the  stress  intensity  factors.  Aberson  and  Anderson  [8] 
used  a  special  crack-tip  singularity  element  to  compute  the  stress  intensity 
factors  for  a  crack  in  a  nonsymmetrical  aft  lug  of  an  engine  pylon.  Pian 
et  al  [9]  used  the  hybrid  finite  element  method  to  compute  the  K-values  for 
cracks  oriented  in  various  angles  from  the  axial  direction  of  straight  lugs. 
Impel lizzeri  and  Rich  [10] modif ied  the  exact  weight  function  derived  by 
Bueckner  [ll],  for  an  edge  crack  in  a  semi-infinite  plate,  to  include  a  scries 
of  geometry  correction  factors.  Then  they  computed  the  K-values  using  the 
weight  function  method.  Except  for  Reference  [8], all  of  these  works  made  the 
assumption  that  the  assumed  or  computed  pin-bearing  pressure  distribution  for 
an  uncracked  case  remains  unchanged  even  after  the  crack  has  initiated  and 
propagated.  Based  on  the  parametric  study  conducted  in  Kef erence  [9], it  was 
found  that,  for  any  given  crack  length,  the  difference  in  the  stress  inten¬ 
sity  factor  computed  using  the  uniform  and  cosine  pin-bearing  pressure  dis¬ 
tributions  was  as  much  as  30  percent.  Therefore,  it  is  salient  that  the 
correct  representation  of  the  pin-bearing  pressure  distribution  during  the 
crack  growth  process  is  essential  to  the  calculation  of  accurate  stress 
intensity  factors.  An  analysis  procedure  using  a  finite  element  method  with 
inclusion  of  a  crack-tip  singularity  element  for  analyzing  cracks  in  both 
straight  and  tapered  lugs  having  a  neat  fit  between  the  pin  and  the  hole[l2] 
has  shown  that  it  can  accurately  account  for  the  change  of  pin-bearing 
pressure  distribution  with  the  change  in  crack  length,  and  provide  stress 
intensity  factors  which  are  in  excellent  agreement  with  the  available  data. 
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This  approach  was  extended  for  the  development  of  analytical  methodology 
for  analyzing  cracks  in  attachment  lugs  with  and  without  the  presence  of 
residual  stress  resulting  from  the  installation  of  the  interference-fit 
bushing  [13].  The  analytical  procedure  developed  for  computing  the  strcr-s 
intensity  factors  for  cracks  in  attachment  lugs  having  residual  stress?;* 
around  the  hole  prior  to  the  application  of  pin-bearing  loads  consist*  of 
two  major  steps.  First,  the  effective  unflawed  stress  distribution  <rn  the 
prospective  crack  surface  was  obtained  by  superimposing  the  residua',  hoop 
stresses  due  to  the  installation  of  an  interference-fit  bushing  on  the 
applied  tangential  stresses  obtained  due  to  the  application  of  a  n'in  loading. 
Second,  a  through  crack  was  introduced  in  this  stress  field  bv  removing  the 
tractions  on  the  crack  faces  and  computing  the  corresponding  stvvss  intensity 
factor  using  the  weight  function  approach. 

This  report  describes  the  analytical  methods  and  procedures  for  obtain¬ 
ing  the  stress  intensity  factors  and  for  predicting  the  fatigue  crack  growth 
life  for  cracks  at  attachment  lugs  with  and  without  the  presence  of  residual 
stresses.  Two  types  of  attachment  lugs  are  considered  in  the  analysis.  They 
are:  (1)  straight-shank  male  attachment  lugs  and  (2)  tapered  male  attachment 
lugs  as  shown  in  Figure  2-1.  Three  different  outer-fo-inner  radius  ratios, 
Ko/R.  =  1.10,  2.25  and  3.0,  are  considered  in  the  analysis.  The  straight 
attachment  lugs  are  subjected  to  axial  pin  loading  only,  while  the  tapered 
attachment  lugs  are  subjected  to  axial,  off-axis  and  transverse  pin  loadings. 
Types  of  cracks  considered  include  single  t hrough-the-t hickness  crack  and 
single  corner  crack  as  depicted  in  Figure  2-2. 

For  through-the-thickness  cracks,  the  methods  developed  and  described 
in  this  report  include:  (1)  the  compounding  solution  method  which  involves 
the  superposition  of  known  solutions  for  idealized  cracked  geometries; 

(2)  the  two-dimensional  cracked  finite  element  method  which  is  capable  of 
characterizing  the  crack-tip  stress  singularity  internally;  and  (3)  the 
weight  function  (also  known  as  Creen's  function)  method  which  computes  the 
stress  intensity  factors  from  the  knowledge  of  the  unflawed  stress  distribu¬ 
tion  and  the  superposition  technique.  For  corner  cracks,  the  methods  devel¬ 
oped  and  described  include:  (1)  the  one  parameter  compounding  method  which 
assumes  a  constant  crack  shape  and  estimates  the  stress  intensity  factor  at 
only  one  location  ou  the  erack  front  (at  lug  surface);  (?)  a  two-dimens ional 
approach  which  estimates  the  stress  intensity  factors  along  the  crack  front 
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(a)  STRAIGHT  ATTACHMENT  LUG 


(b)  TAPERED  ATTACHMENT  LUG 


KUure  2-1 .  types  ot  Attachment  t 

Considered  in  the  Analyst* 


of  a  corner  crack  by  modifying  the  Green's  function  for  through-the-thickness 
crack  solutions  with  appropriate  hole  curvaturev  flaw  shape,  front  free 
surface  and  back  surface  correction  factors;  (3)  the  rigorous  three- 
dimensional  cracked  finite  element  method,  which  was  developed  using  a 
hybrid  displacement  finite  element  procedure. 

After  the  stress  intensity  factor  is  determined  as  a  function  of  crack 
length  for  a  particular  lug  and  loading  condition,  its  range  (AK)  can  be 
computed  for  the  Nth  cycle  in  a  given  load  spectrum.  Then  a  numerical  inte¬ 
gration  of  da/dN  versus  AK  relation  can  be  carried  out  to  predict  the  crack 
growth  characteristics.  An  automated  computer  program  was  developed  utiliz¬ 
ing  the  state-of-the-art  fracture  mechanics  methodologies  for  the  prediction 
of  fatigue  crack  growth  history  under  fatigue  loading.  The  program  is  cap¬ 
able  of  predicting  the  crack  growth  behavior  of  single  corner  cracks  and 
single  through-the-thickness  cracks  in  attachment  lugs  using  block-by-block 
integration  technique.  In  the  case  of  corner  crack  problems,  a  transitional 
crack  growth  criterion  from  part-through  crack  to  through-the-thickness 
crack  was  also  developed  and  incorporated.  Several  load-interaction  models 
were  also  included  in  the  crack  growth  analysis  program  for  spectrum  load¬ 
ings. 

The  prediction  methods  presented  in  this  volume  of  the  report  are  to  be 
evaluated  against  fatigue  crack  growth  test  results.  The  test  program  con¬ 
sists  of  two  parts.  Group  l  tests  consist  of  straight-shank  male  attachment 
lugs.  Both  constant  amplitude  and  block  spectrum  loadings  are  considered  in 
this  group.  Variational  tests,  such  as  different  thickness  and  with  residual 
stresses  due  to  interference-fit  bushings,  are  also  included  in  Group  1  test¬ 
ing.  Group  It  tests,  performed  later,  encompass  several  additional  lug 
configurations  (tapered,  dogbone,  clevis  and  a  simulated  wing-pylon  attach  lug) 
and  include  both  axial  and  off-axis  loadings. 

The  stress  intensity  factor  formulations  prepared  for  Croup  I  straight- 
shank  male  attachment  lugs  are  described  in  Section  III.  A  draft  of  Section 
III  was  completed  prior  to  Gro^  I  testing.  After  Graf>  I  testing,  these 
formulations  are  evaluated.  Based  on  this  evaluation,  a  stress  intensity 
factor  estimation  methodology  is  selected  for  Group  II  testing.  Section  IV 
describes  such  a  selected  methodology  for  tapered  attachment  lugs.  The 
crack  growth  prediction  method  which  employs  these  stress  intensity  factor 
solutions  is  presented  in  Section  V. 
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SECTION  III 


STRESS  INTENSITY  FACTORS  FOR  STRAIGHT  ATTACHMENT  LUGS 


In  this  section,  the  analytical  methodologies  required  for  the  develop¬ 
ment  of  stress  intensity  factors  for  cracks  in  simple  straight-shank  male 
attachment  lugs  are  described.  The  stress  intensity  factor  solutions  are 
obtained  for  single  through-the-thickness  cracks  and  single  quarter- 
elliptical  corner  cracks  in  attachment  lugs.  As  mentioned  previously,  the 
following  methods  are  used  for  through-the-thickness  and  corner  cracks. 

For  through-the-thickness  crack  problems: 

(1)  Simple  compounding  method 

(2)  Two-dimensional  cracked  finite  element  procedure 

(3)  Weight  function  or  Green's  function  method 

For  corner  crack  problems: 

(1)  Simple  one  parameter  compounding  method 

(2)  Two-dimensional  Green's  function  solution  modification 

(3)  Rigorous  three-dimensional  cracked  finite  element  procedure 

This  section  also  describes  the  method  for  developing  the  stress  in¬ 
tensity  factor  for  problems  with  residual  stresses  (due  to  the  installation 
of  interference-fit  bushings).  A  stress  analysis  of  unflawed  straight 
attachment  lugs  was  also  performed  and  is  presented  here.  Note  that  all 
the  specimens  tested  in  the  Group  I  test  program  are  straight  attachment 
lugs  and  the  methodologies  described  in  this  section  are  evaluated  using 
Croup  I  test  results. 

L. _ STRESSES  IN  AN  UN FLAWED  STRAIGHT  ATTACHMENT  LUC 

A  primary  item  of  interest  in  the  analysis  of  cracked  attachment  lugs 
is  the  stress  distribution  in  the  uncracked  state.  This  is  important  for  at 
least  three  reasons.  The  location  of  the  peak  tangential  stress  determines 
the  location  of  the  most  critical  crack.  The  stress  intensity  factor  for  a 
very  small  crack  is  proportional  to  the  stress  concentration  factor  value. 
Finally,  the  weight  function  method  of  estimating  stress  intensity  factors 
requires  the  stresses  on  the  respective  crack  surface  in  the  uncracked  lug. 
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The  conventional  displacement  finite  element  method  was  used  to  obtain, 

for  an  uncracked  straight  lug,  (1)  the  tangential  stresses  along  the  edge  of 

the  hole,  (2)  the  stress  distribution  on  the  prospective  crack  surface,  and 

(3)  the  pin-bearing  pressure  distribution  along  the  contact  surface  between 

the  pin  and  the  lug.  The  rigidity  ratio,  E  .  /E,  ,  is  assumed  to  be  3.0  in 

pin  lug 

the  present  analysis. 

Figure  3-1  shows  a  typical  model  used  in  the  analysis.  Due  to  synmetry, 
only  the  upper  half  of  the  lug  was  modeled.  The  lug  and  the  pin  were  rep¬ 
resented  by  a  set  of  constant-strain  triangular  and  quadrilateral  elements.- 
To  load  the  model,  a  concentrated  force  was  applied  at  the  center  of  the  pin 
and  reacted  at  the  other  end  of  the  lug.  Spring  elements  (S)  were  used  to 
connect  the  pin  and  lug  at  each  pair  of  nodes  having  identical  nodal  coor¬ 
dinates  all  around  the  periphery.  The  area  of  contact  was  determined  itera¬ 
tively  by  assigning  a  very  high  stiffness  to  spring  elements  which  were  in 
compression  and  a  very  low  stiffness  (essentially  zero)  to  spring  elements 
which  were  in  tension.  A  neat-fit  and  no  friction  assumptions  were  made  for 
the  pin-to-lug  assembly. 

The  calculated  tangential  stresses  (normalized  by  the  average  bearing 
stress)  along  the  edge  cf  the  hole  for  a  pin  loading  applied  in  the  0°  direction 
are  shown  in  Figure  3-2  for  RQ/Rl  ratios  of  1.50,  2.25  and  3.00.  As  antici¬ 
pated,  the  maximum  tangential  stresses  are  located  at  about  90°  away  from 
the  axis  of  the  lug,  and  the  minimum  stresses  (small  amount  in  compression ) 
are  located  at  the  180°  location. 

The  tangential  stresses  for  pin  loading  applied  in  axial  compression 
(0  -  180  )  are  shown  in  Figure  3-3.  As  can  be  seen  from  this  figure,  the 
resulting  stresses  at  0  =  >90  are  small  tensile  stresses.  Thus,  reversing 
the  loading  direction  does  not  cause  a  stress  reversal  at  the  critical  90° 
location.  Furthermore,  the  tensile  stresses  produced  are  relatively  small. 

The  unflawed  elastic  stress  distributions  on  the  prospective  crack 

plane  for  axial  tension  applied  in  the  0°  direction  are  presented  in  Figure 

3-4.  These  values  are  also  tabulated  in  Table  3-1.  The  tangential  stresses 

along  the  x-axls  shown  in  the  figure  are  normalized  by  the  average  bearing 

stress.  It  is  clear  that  the  gradient  of  the  stress  distribution  close  to 

the  edge  of  the  hole  is  very  steep,  especially  for  a  lug  having  a  smaller 

R^/R^  ratio.  The  corresponding  elastic  stress  concentration  factors  at  the 

edge  of  the  hole,<?  /o  ,  are  plotted  in  Figure  3-5.  The  equation  of  the 
ttu^x  b  r 
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logarithmic  straight  line  shown  in  this  figure  fits  the  computed  values  with¬ 
in  0.5  percent,  as  shown  in  the  accompanying  table. 

The  computed  pin-bearing  pressure  distributions  along  the  contact  sur¬ 
face  are  shown  in  Figure  3-6.  The  figure  shows  that  the  shape  of  the  dis¬ 
tribution  is  close  to  uniform  at  the  central  portion  of  the  contact  surface. 

As  the  R  /R.  ratio  increases,  the  pressure  decreases  near  0  =  +  90°  but 

o  i  o  — 

increases  near  0=0. 

2.  STRESS  INTENSITY  FACTORS  FOR  THROUGH-THE-THICKNESS  CRACKS 

The  development  of  stress  intensity  factor  solutions  for  single  through- 
the-thickness  cracks  in  attachment  lugs  by  various  methods  are  described  in 
this  subsection.  Also  included  in  this  subsection  are  the  methods  for  analyz¬ 
ing  attachment  lugs  with  interference-fit  bushings  and  elasto-plastic  analy¬ 
sis  for  analyzing  lugs  which  are  loaded  above  yield. 

2 . 1  COMPOUNDING  METHOD  APPROXIMATIONS 

Stress  intensity  factors  for  an  axially-loaded  straight  lug  with  a 
single  through-the-thickness  crack  can  be  estimated  by  the  compounding  of 
known  K  solutions.  The  compounding  method  has  been  used  for  crack  growth 
analysis  of  complex  structure  [14]. 

As  Figure  3-7  shows,  the  actual  lug  configuration  is  approximated  as  an 
infinite  strip  of  width  2  Ro  containing  a  central  hole  of  radius  R^.  A 
single  through-the-thickness  crack  is  assumed  to  be  present  at  9  =  90°  to 
the  axis  of  the  lug.  The  length  of  the  crack  is  a^,  measured  from  the  edge 
of  the  hole.  The  total  pin  load  P  is  assumed  to  result  from  a  uniform 
distribution  of  radial  pressure  from  0  =  -67.5°  to  67.5°,  and  to  be  reacted 
by  a  remote  uniform  stress  (j q  =  P/(2  K^B),  where  B  is  the  lug  thickness. 

The  stress  intensity  factor  for  the  configuration  shown  in  Figure  3-7 
can  be  approximated  by  combining  a  number  of  known  K  solutions.  The  basic 
solutions  are  for  cracks  emanating  from  holes  in  infinite  plates.  Multi¬ 
plicative  correction  factors  are  used  to  modify  these  infinite  plate 
solutions  to  account  for  the  effects  of  finite  width. 
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2.1.1  Cracks  at  Holes  In  Infinite  Plates  -  The  loading  shown  in  Figure  3-7 
consists  of  pin  loading  in  one  direction  and  remote  loading  in  the  other. 

Tweed  and  Rooke  [15]  have  obtained  the  solution  for  a  crack  emanating  from 
a  circular  hole  in  an  infinite  plate  subjected  to  remote  loading.  This  is 
Case  31  in  the  compendium  given  in  Appendix  A.  The  formula  for  K31  fits 
Tweed  and  Rooke' s  solution  within  0.5  percent. 

The  general  solution  for  the  same  geometry  subjected  to  arbitrary  radial 
loading  at  the  hole  is  given  in  Reference  [16].  Case  35  in  Appendix  A  gives  an 
equation  which  closely  fits  the  special  case  of  uniform  radial  pressure  ap¬ 
plied  from  8  =  -67.5°  to  0=  67.5°,  and  from  (180°  -  9)  =  -67.5°  to 
(180°  -  8)  =  67.5°. 


2.1.2  Symmetrical  Line  Cracks  in  Strips  -  The  solution  for  a  symmetrical 

line  crack  of  length  2a^  in  an  infinite  strip  of  width  2b  under  uniform 

tension  stress  <jo  is  well  known.  It  is  given  in  Appendix  A  as  Case  21. 

The  function  $,,(X  )  is  called  the  width  correction  factor,  because  it  is 
21  o 

applied  multiplicatively  to  the  infinite-plate  K  solution  to  give  the 
correct  solution  for  the  strip. 

Case  22  in  Appendix  A  gives  an  approximate  solution  for  the  symmetrically- 
cracked  strip  subjected  to  concentrated  splitting  forces  P  applied  to  both 
surfaces  of  the  crack  at  the  centerline.  This  approximation  is  derived  by  the 
similarity  method  from  three  known  K  solutions,  as  followsi 

Figure  3-8  shows  the  four  similar  configurations.  In  all  four  cases  the 
crack  lengths  are  the  same  and  the  shear  stresses  vanish  along  the  lines 
x  -  _»b.  The  unknown  solution  is  Kj2*  solutions  for  the  other  three  cases 
shown  are  known  in  closed  form.  is  the  well-known  solution  for  periodic 

col  linear  cracks  subjected  to  uniform  loading  [17]: 
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The  solution  for  periodic  col  linear  cracks  subjected  to  central 

splitting  forces  on  the  crack  surfaces  is  given  by  Tada,  et  al  [18]. 
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Notice  in  Figure  3-8  that  Cases  21  and  22  are  identical  geometries,  as 
are  Cases  19  and  20.  Furthermore,  Cases  20  and  22  are  identical  in  the  ap¬ 
plied  loadings,  as  are  Cases  19  and  21. 

Thus,  the  correction  factor  approach  can  be  applied  in  two  alternative 
ways  to  estimate  The  ratio  (Kj^/K^),  a  correction  factor  to  account 

for  the  difference  between  a  finite-width  strip  and  a  periodic  crack  array, 
can  be  multiplied  times  K^q,  the  solution  for  the  periodic  crack  array  with 
the  Case  22  loading  condition.  Alternatively,  the  ratio  (K^q/K^),  a  correc“ 
tion  factor  to  account  for  the  difference  between  remote  loading  and  point 
loading  on  the  crack  line,  can  be  multiplied  times  the  solution  for 

remote  loading  of  the  Case  22  geometry.  Either  way,  the  resulting  equation 
for  1^2  is  the  same: 


K21  *20 
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2.1.3  Effect  of  Eccentricity  -  The  crack  in  the  lug  is  eccentrically  located, 
so  eccentricity  must  be  considered.  Isida  [19]  has  obtained  the  stress  in¬ 
tensity  solution  for  an  unsymmetrical  crack  in  a  strip  under  uniform  tension. 
Equations  are  given  in  Appendix  A,  Case  26,  which  fit  Isida's  numerical 

results  within  about  3  percent  (within  1  percent  for  a/b  <  0.6).  There  are 

(N)  n 

two  stress  intensity  formulas,  1C.  for  the  crack  tip  nearest  the  edge  and 
(F)  26 

for  the  crack  tip  farthest  from  the  edge.  For  the  symmetrical  crack 
both  formulas  reduce  to  the  equation  for  Kjj* 

Case  2 7  in  Appendix  A  gives  an  approximate  solution  for  the  unsymmetri¬ 
cal  crack  subjected  to  splitting  forces  on  the  crack  surfaces  at  the  strip 
centerline.  This  solution  is  obtained  by  the  Similarity  Method  from  the  five 
known  K  solutions  shot.  .  in  Figure  3-9.  (The  solutions  for  Cases  12  and  13 
are  well  known  and  are  given  in  deference  [ 18].)  There  are  two  eccentricities 
in  Case  27:  the  crack  is  not  centered  on  the  strip  centerline,  and  the  load 
is  not  centered  on  the  crack  centerline.  Remove  both  eccentricities  and 
Case  22  is  obtained. 
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The  crack  is  not  centered  on  the  strip  in  Case  26.  Remove  this  eccentri¬ 
city  and  Case  21  is  obtained.  Therefore,  the  correction  factor  to  account 

for  eccentricity  between  the  crack  and  the  strip  is  the  ratio  ^6^2 1  * 

The  load  is  not  centered  on  the  cr/ick  in  Case  13.  Remove  this  eccentri¬ 
city  and  Case  12  is  obtained.  Therefore,  the  correction  factor  to  account 

for  eccentricity  between  the  loading  po.  at  and  the  crack  is  the  ratio  K^/K^ 

When  both  correction  factors  are  applied  to  the  approximate  formula 

for  K^,  8*ven  Appendix  A,  is  obtained: 


,(N)  „(N) 
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(4) 


2.1.4  Cracks  at  Holes  in  Strips  -  The  stress  intensity  solutions  for  cracks 
at  holes  in  infinite  plates  and  for  line  cracks  in  strips  can  be  combined  to 
obtain  approximate  K  formulas  for  cracks  at  holes  in  strips.  The  resulting 
formulas  can  be  superimposed  to  estimate  K  for  a  pin-loaded  lug  with  a 
crack. 

Case  41  in  Appendix  A  gives  an  approximate  solution  for  a  uniformly 
loaded  strip  with  a  crack  emanating  from  a  central  hole.  Referring  to 
Figure  3-10,  one  would  expect  the  solution  for  to  be  approximated  by 

the  product  *26^11* 

This  product,  however,  does  not  account  for  the  interaction  between  the 
hole  and  the  edges  of  the  strip.  The  gross  area  stress  concentration  factor 
given  by  Peterson  [  20 ]  for  a  center  hole  in  a  strip  can  be  approximated  by 
the  following  formula: 


kt(r/b)  -  3  sec  (~)  -  3  4^  (r/b)  (S) 

As  the  crack  length  a ^  approaches  zero,  has  a  known  asymptotic 
solution  in  terms  of  the  stress  concentration  factor,  i.e. 

Si  2 

Urn  (K4l)  -  kt  (r/b)  -4^  (r/b)  (6) 

aj/r  .0 
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To  satisfy  Equation(6),  a  hole-strip  width  correction  factor  <J>21 
must  be  included,  so  that  the  approximate  solution  for  Case  41  is 


St  K26 
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(8) 


Case  44  in  Appendix  A  is  an  approximate  solution  for  a  strip  with  a 
crack  at  a  central  hole,  loaded  synroetrically  at  the  hole.  This  solution  is 
derived  by  the  same  approach  used  to  derive  K41‘  In  this  case,  however, 
the  gross  area  stress  concentration  factor  for  the  uncracked  case  is  not 
known.  Figure  3-5,  obtained  by  finite  element  analysis,  could  have  been  used 
here  to  estimate  the  required  k£,  but  to  do  so  would  have  compromised  the 
independence  of  the  solutions  generated  by  the  compounding  and  finite  element 
methods.  Based  on  engineering  judgment,  therefore, <t>2^  (r/b)  is  appended  to 
product  solution  obtained  from  Figure  3-11  to  approximate  K^,  just  as  was 
done  for  Case  41.  Thus,  is  approximated  by  the  equation 


44 


K35  K27 
K13 


^21 


(9) 


2.1.5  Cracked  Lug  -  The  final  step  in  the  derivation  of  a 
method  stress  intensity  factor  for  the  pin-loaded  straight 
superposition.  As  demonstrated  in  Figure  3-12,  the  result 
Caes  41  and  44  is  identical  to  the  result  of  superimposing 
nation  in  Figure  3-7  with  itself  (oriented  upside-down  and 


compounding 
lug  is  a  simple 
of  superimposing 
the  lug  approxi- 
backwards).  Thus 


Kuu=  •  I 


(10) 


Equation  (10)  was  used  to  compute  the  stress  intensity  factors  for 
R  /R.  ratios  of  1.50,  2.25  and  3.0,  and  the  results  are  tabulated  in 

O  l 

Table  3-2. 


2.2  TWO-DIMENSIONAL  CRACKED  FINITE  EhSMSNT  METHOD 

Because  of  the  ease  with  which  the  f in  it e-element  method  handies  com¬ 
plex  geometries  and  boundary  conditions,  chis  method  has  been  used  exten¬ 
sively  to  study  fracture  in  complex  structures.  Two  special  crack-tip 
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for  Approximating 


29 


6.287  1.659  0.851  9.431  3.733  2.553 


singularity  elements,  usually  referred  to  as  cracked  elements,  developed  at 
Lockheeri-Georgia  [8]  have  brought  the  power  and  flexibility  of  the  finite- 
element  method  to  bear  much  more  effectively  on  fracture  mechanics  problems. 
This  method  has  been  used  to  analyze  a  crack  in  an  unsymmetric  lug  of  a 
C-5  engine  pylon  [8].  It  has  been  modified  slightly  to  accurately  account 
for  the  change  in  pin-bearing  pressure  distribution  for  the  change  in  crack 
length.  The  following  discussion  describes  this  method  as  it  was  used  in 
the  calculation  of  stress  intensity  factors  for  cracks  at  attachment  lugs. 

Figure  3-13  shows  a  typical  finite  element  model  used  for  a  single  crack 
emanating  from  a  straight  attachment  lug.  A  10-node  high-order  singularity 
element  [8  ]  was  used  at  the  crack  tip  region  for  computing  the  stress  inten¬ 
sity  factors.  The  remainder  of  the  finite  element  modeling  details  was 
identical  to  those  described  in  the  stress  analysis  section  for  the  uncracked 
lug,  including  the  iterative  determination  of  the  contact  area.  The  solution 

was  obtained  for  several  outer-to-inner  radius  ratios,  R  /R, ,  and  for  a 

o  i 

succession  of  cracks  having  normalized  lengths,  c/(Rq  -  R^,)  ranging  from  0.1 
to  0.9.  As  in  the  case  of  stress  analysis,  the  rigidity  ratio,  , 

of  3.0  is  assumed  in  the  fracture  analysis. 

The  pin-bearing  pressure  distributions  obtained  for  single  cracks  ema¬ 
nating  from  attachment  lugs  loaded  by  neat-fit  pins  are  presented  in  Figures 
3-14  through  3-16.  These  figures  show  the  effect  of  crack  length  on  the  pin¬ 
bearing  pressure  distributions  for  a  straight  lug  having  Rq/R^  ratios  of 
1.5,  2.25  and  3.00,  respectively.  When  there  is  no  crack,  the  distribution 
is  close  to  uniform  at  the  central  portion  of  the  contact  surface,  unlike 
the  cosine  distribution  commonly  assumed  in  literature.  There  is  no  contact 
at  the  expected  crack  location,  0  =  -90°.  However,  as  soon  as  a  crack 
appears,  the  contact  spreads  to  one  side  ("upper  lip")  of  the  crack  mouth. 

As  the  crack  opens,  the  other  side  ("lower  lip")  of  the  crack  mouth  moves 
away  from  the  pin  resulting  in  no  pin-bearing  pressure  on  the  lower  lip. 
Figures  3-14  through  3-16  show  that  the  longer  the  crack,  the  higher  the 
pin-bearing  pressure  on  the  upper  lip  of  the  crack  mouth,  especially  for  a 

small  R  /R.  ratio.  For  example,  for  a  R  /R,  ratio  of  1.5,  the  normalized 
o  i  o  i 

pressure  at  the  upper  lip  of  the  crack  mouth  increases  from  zero  for 
c/R^  —  0  to  more  than  twice  the  average  pressure  for  c/R^  >  0.3.  The  pres¬ 
sure  distribution  elsewhere  on  the  contact  surface  also  changes  with  the 
crack  length.  For  larger  crack  sizes,  the  pressure  decreases  from  an  initial 
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maximum  value  at  the  crack  mouth  to  a  local  minimum  and  then  gradually  in¬ 
creases  again  and  approaches  another  maximum  before  decreasing  to  zero  at 
the  end  of  the  contact  surface. 

The  effect  of  pin-bearing  pressure  distributions  on  stress  intensity 
factors  is  shown  in  Figure  3-17  for  single  cracks  in  a  straight  lug  having 
a  R^K^  ratio  of  1.5.  In  Figure  3-17,  the  circle  and  triangle  symbols  are 
results  obtained  from  Reference  [9]*,  where  it  is  assumed  tha .  the  pin-bearing 
pressure  distributions  are  cosine  and  uniform  distributions  along  the  180- 
degree  contact  surface,  respectively,  and  that  such  distributions  remain  un¬ 
changed  with  the  crack  length.  The  square  symbol  represents  the  results 
obtained  using  the  present  analysis,  which  has  properly  accounted  for  the 
change  in  pressure  distribution  as  the  result  of  crack  extension.  As  can  be 
seen  from  this  figure,  when  the  crack  is  small  (c/R^  <  0.05),  the  current 
computed  K-value  is  practically  the  same  as  that  obtained  using  a  uniform 
pin-bearing  pressure  distribution.  However,  as  the  crack  length  increases, 
the  current  analysis  gives  a  lower  K-value  than  the  others.  This  is  because 
when  the  crack  length  increases,  the  pressure  near  Q  =  -90°  increases 
markedly  and  exceeds  the  average  pin-bearing  pressure  (see  Figure  3-14). 

This  high  pin-bearing  pressure,  when  applied  near  9  -  -90°  in  the  direction 
almost  parallel  to  the  crack  orientation,  tend  to  close  the  crack  surfaces, 
hence  reducing  the  stress  intensity  factor  as  discussed  by  Brussat  [16]. 

The  computed  normalized  opening  mode  stress  intensity  factors  using  a 

steel  pin  and  an  aluminum  lug  model  are  presented  in  Figure  3-18  as  a  function 

of  normalized  crack  length  (c/K^)  for  single  cracks  emanating  from  the  hole 

wall  of  straight  attachment  lugs  with  ^/R^  ratios  ranging  from  1.5  to  3.0. 

In  all  cases,  the  computed  sliding-mode  stress  intensity  factors,  are 

much  smaller  than  those  of  the  opening  mode,  K^ ,  so  they  are  not  presented 

in  the  figures.  It  should  be  noted  that  the  K-vaiues  were  normalized  in 

terms  of  the  average  bearing  stress,  <7^,  instead  of  far-field  gross  section 

stress,  <j  .  To  convert  these  normalized  factors  in  terms  of  Q  ,  one  can 
o  o 

simply  multiply  these  normalized  factors  by  the  corresponding  ratio  of 
Ko/k . .  For  convenience,  the  computed  K-values  normalized  in  terms  of  the 
average  far-ftcld  stress  are  tabulated  in  Table  3-3  as  a  function  of 

*The  results  are  mislabeled  In  Reference  [9];  uniform  distribution  results  are 
labeled  cosine  distribution  and  vice  versa. 
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normalized  crack  length  (normalized  by  the  net  section).  Note  that 
at  the  edge  of  the  hoie  (c/R^  ■  0),  the  normalized  stress  intensity  factor 
was  obtained  by  multiplying  the  stress  concentration  factor  determined 
from  unflawed  stress  analysis  by  1.12,  which  was  derived  by  Gross  et  al  [21  ] 
for  a  straight  edge  crack  in  a  finite-plate  specimen  loaded  in  tension.  Un¬ 
flawed  elastic  distributions  on  the  prospective  crack  plane  and  the  corre¬ 
sponding  stress  concentration  factors  at  the  edge  of  the  hole  have  been  given 
in  Figures  3-4  and  3-5,  respectively. 


Studies  were  also  made  to  investigate  the  effects  of  a  single  crack 
versus  a  symmetrical  double  crack  and  relative  rigidity  of  the  pin  and  the 
lug  (Ep^/Ej^).  A  typical  comparison  of  K-values  obtained  for  a  single 
crack  and  a  symmetrical  double  crack  in  a  straight  aluminum  lug  loaded  by  a 
steel  pin  (E  3  3.0)  is  shown  in  Figure  3-19.  As  anticipated,  the 

computed  K-values  for  a  double  crack  are  higher  than  those  of  a  single  crack, 
and  the  difference  increases  with  the  crack  length.  It  increases  from  less 
than  1Z  for  c/R^  <  0.2  to  about  37Z  for  c/H^  =  0.9.  To  study  the  effect 
of  the  relative  rigidity  of  the  pin  and  the  lug  on  the  stress  intensity 
factors,  the  computed  K-values  for  a  single  crack  emanating  from  the  hole 
wall  of  a  steel  lug  loaded  by  a  steel  pin  (B  s  1-0)  15  also  included 

in  Figure  3-19.  As  presented  in  this  figure,  the  K-values  computed  for  a 
combination  of  a  steel  pin  and  steel  lug  is  slightly  higher  than  that  of  the 
steel  pin  and  the  aluminum  lug.  A  comparison  of  stress  distributions  along  the 
0  *  -90°  radial  line  for  steel  lug-steel  pin  and  aluminum  lug-steel  pin 
combinations  is  presented  in  Figure  3-20.  These  stresses  can  be  used  in 
conjunction  with  Creen's  function  Co  calculate  the  stress  intensity  factors 
appropriately. 


£.3  THE  WKIGH1  FUNCTION  METHOD 

The  linear  superposition  method  has  been  used  frequently  to  obtain  the 
stress  intensity  factors  for  various  types  of  crack  problems.  The  principle 
of  superposition  of  linear  elasticity  implies  chat,  lor  the  purpose  of  cal¬ 
culating  stress  intensity  factors,  loading  the  crack  faces  with  or(x)  is 
equivalent  to  loading  the  cracked  body  with  loads  which  produce  (J(x)  on 
the  prospective  crack  faces  in  the  absence  of  a  crack. 

Figure  3-21  shows  the  scheme  of  the  linear  superposition  method.  The 
stress  intensity  factor  of  problem  3-?la  is  equivalent  to  the  sum  of  that  of 
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lln f  lawed  Stress  Distributions  in  Straight  Attachment  f. 

lor  Various  E  /E,  Ratios  (R  /R.  =  2.25) 
pin  lug  o  i 


problems  3-21b  and  3-21c.  Since  problem  3-21b  is  crack  free,  the  stress- 
intensity  factor  of  problem  3-21a  is  equivalent  to  that  of  problem  3-21c. 

By  idealizing  the  stress  in  problem  3-21c  as  N  discrete  loads,  P  ....,P  , 
then  the  stress-intensity  factor,  for  a  given  crack  length  c,  can^e  com¬ 
puted  from  the  following  equation 

N  N 


K(c)  -  ^  ^  ki(xi,c)  Pt  (*t) 


(11) 

th 


where  k^x^c)  is  the  normalized  stress-intensity  factor  due  to  the  i 
load,  Pj,  applied  at  location  x^  For  arbitrary  distributed  stress,  a  (x) 
instead  of  discrete  forces,  P£,  Equation  (11)  becomes 


K(c)  » 


c 

J  k(x,c)  . 


<7(x)  dx 


(12) 


In  Equation  (12),  k(x,c)  is  the  weight  function  (or  Green’s  function). 
Bueckner  [22]  and  Rice  [23]  defined  the  weight  function  as 


k(x,c) 


H  du(x,c) 

Me)  *  a~ 


(13) 


for  a  symmetrical  load  system  on  a  linearly  elastic  body  containing  a  crack 
of  length  c.^ln  the  above  equation,  H  is  an  appropriate  elastic  modulus: 
it  is  E/(i_p  )for  plane  strain  and  E  for  generalized  plane  stress.  K(c) 
is  the  known  stress  intensity  factor  and  u(x,c)  is  the  y-component  of  the 
crack  surface  displacement  at  x.  The  weight  function  was  shown  [22,23]  to 
be  unique  for  a  given  structural  geometry  and  crack  size  regardless  of  the 
loading  condition.  Therefore,  it  can  be  developed  for  one  load  condition 
and  then  utilized  to  determine  the  stress-intensity  factor  for  any  other 
load  condition. 

The  closed  form  expressions  for  the  weight  function  for  edge  cracks 
[U.24],  center  cracks  [25]  and  coll  inear  cracks  [26]  in  a  wide  panel  are 
available.  However,  the  closed  form  weight  function  for  cracks  emanating 
from  a  hole  is  not  available.  Therefore,  the  weight  function  for  a  straight 
center  crack  has  sometimes  been  used  to  estimate  the  stress-intensity  factor 
£or  radial  cracks  emanating  from  a  circular  hole  [27-29  ],  located  in  the 
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geometric  center  of  a  long  plate.  Impellizzeri  and  Rich  [to]  modified  the 
weight  function  derived  by  Bueckner  [ll]  for  an  edge  crack  in  a  semi-infinite 
plate  to  include  the  geometric  correction  factors  for  estimating  the  stress- 
intensity  factor  of  a  crack  in  an  attachment  lug.  In  order  to  establish  the 
appropriate  values  of  the  derivative  du/dc,  to  be  used  in  Equation  (13), 

Grandt  [30 ]  supplemented  Bowie's  solutions  [31 ]  for  a  radially-cracked  hole 
with  a  finite  element  analysis  to  obtain  the  corresponding  crack-mouth  open¬ 
ing  displacements.  Hsu  and  Rudd  [32]  developed  thd  Green's  function  from 
the  stress-intensity  factors  computed  using  the  high-order  singularity  ele¬ 
ment  for  a  double-radial  crack  emanating  from  an  open  hole  and  subjected  to 
a  pair  of  concentrated  loads  on  and  perpendicular  to  the  crack  surface. 

This  approach  is  capable  of  accounting  for  the  effect  of  hole  curvature, 
finite  width  and  the  profile  of  the  lug  head.  Therefore,  it  was  used  to 
develop  the  weight  function  for  through-the-thickness  cracks  emanating  from 
the  hole  of  attachment  lugs. 

1/2 

By  defining  the  weight  function  G  =  k(c/7T)  and  i;  =  (x-R^)/c  and  sub¬ 
stituting  them  into  Equation  (12),  one  obtains 

K(c)  -  <7o  ViFc  f  G(  c,  £)  a  (£)  d£  (14) 

o 

where  nr  is  the  uniform  far-field  stress  and  <J  =  <j/(T  is  the  normalized 
o  o 

unflawed  stress  distribution  on  the  prospective  crack  surface. 

The  weight  function  G,  for  a  single  radial  crack  emanating  from  a 

circular  hole  of  the  lug  and  subjected  to  a  pair  of  concentrated  forces  on 

the  crack  surfaces,  as  shown  in  Figure  3-22a,  is  obtained  from  the  computed 

stress-intensity  factor  using  2-D  cracked  finite  element  analysis  for  various 

R  / R,  ratios,  crack  lengths  c/(R  -R.),  and  (x-R,)/c  ratios  as  follows: 
o  i  o  l  A 

R 

G  ((R  -r"V  c  ’  ir^  =  k(x>c)  VaTtf  (15) 

Finite  element  methodology  becomes  difficult  to  employ  when  the  concen¬ 
trated  forces  are  applied  close  to  the  crack  tip,  say  (x-R^)/c  >  0.9. 

In  this  range  the  corresponding  weight  functions  were  obtained  using 

the  edge  crack  model  as  shown  in  Figure  3-22b.  The  weight  function,  G^ ,  for 

such  edge  cracks  is  available  in  Reference  [18],  and  can  be  written  as 


(a)  FOR  P  APPLIED  AWAY  FROM  (b)  FOR  P  APPLIED  CLOSE  TO 

CRACK  TIP  CRACK  TIP 

Figure  3-22,  Model  for  Generating  th  ’  Weight  Functions 


+ 


3.52(l-(x-R.)/c) 
(1  -  c/2b)3/2 


4.35  -  5.28  (x-R^/c 
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r  1.30-0.30  ((x-R.)/c)3/2 

1 

T  r  1 

) 

+  0.83-1.76  (x-R^/c 

1— (1— (x— R  )/c)  l 

J  (16) 

L  \1  -  ((x-R.)/c) 

L 

) 

In  the  conventional  finite  element  method,  the  external  force  can  only 
be  applied  at  nodal  points.  When  the  crack  length  is  small,  it  becomes 
cumbersome  to  refine  the  model  such  that  there  will  be  enough  nodes  along 
the  crack  faces  for  the  purpose  of  computing  the  K  and  G  values.  There¬ 
fore,  an  alternate  approach  is  used.  For  each  crack  length  c,  the  K  and 
G  values  were  calculated  at  each  available  nodal  point  on  the  crack  face, 
say  (Xj-R^l/c,  using  a  finite  element  model  for  the  configuration  shown  in 

Figure  3-22a  and  Equation  (15).  The  weight  function  G^  for  an  edge  crack 
in  a  finite  width  strip  (Figure  3-22b)  was  then  calculated  at  the  same  loca¬ 
tions,  (xi-Ri>/c,  using  Equation  (16),  and  the  ratio  was  obtained. 

r  ((x.-R.)/c)  =  G  ((xj-R.)/c)/G1  ((x^R^/c)  (17) 


Assuming  that  r  approaches  1  as  (x  .-R  .)/c  approaches  1,  a  least 
squares  polynomial  fit  is  obtained. 


M 

r  ((x-R.)/c)  =  £  Cj((x-R.)/c)j_1 


(18) 


The  weight  function  at  any  location  on  tne  entire  crack  surface  can 
then  be  computed  using  Equations  (16)  and  (18)  as 

G  ((x-Rt)/c)  «  r  ((x-Rt)/c)  Cj  ((x-R^/c)  (19) 

The  results  obtained  using  the  aforementioned  procedure  are  shown  in 
Figures  3-23  through  3-25  for  R^/R^  ratios  of  1.50,  2.25,  and  3.00, 
respectively.  The  symbols  shown  in  these  three  figures  are  the  discrete 
values  obtained  using  the  finite  element  method  and  the  curves  are  obtained 
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using  Equation  (19).  As  seen  from  these  figures,  this  alternate  approach 
gives  an  excellent  estimation  of  G-values  at  any  desired  location. 

With  a  knowledge  of  the  weight  functions,  G,  and  the  unflawed  stress, 
g,  on  the  prospective  crack  surface  with  the  crack  absent,  one  can  numeri¬ 
cally  integrate  in  Equation  (14)  to  obtain  the  corresponding  stress-intensity 
factor  for  a  radial  crack  in  an  attachment  lug. 

Integration  of  Equation  (14)  can  be  carried  out  accurately  and  effici¬ 
ently  by  the  Gaussian  quadrature  formula.  In  order  to  perform  the  Gaussian 
integration,  the  integral  in  Equation  (14)  needs  to  be  transferred  into  the 
well-known  Newton-Cotes  quadrature  formula  for  the  range  of  integration  be¬ 
tween  -1  and  +1.  This  can  be  easily  done  by  assuming  that 

2£  -  1,  or  $  =  |  (?>  1),  (20) 


and  the  integral  of  Equation  (14)  becomes 


U)G(c,0d£ 


c 1  n 

/  P(?)d5  =  £  Hi  p(5t) 

-i  i-i 


(21) 


where 


HO  =  |  a  G  (c,  (22) 

The  abscissae  (^)  and  weight  coefficients  (H^)  of  the  Gaussian  quad¬ 
rature  formula,  Equation  (21),  for  n  up  to  10  are  given  in  Reference  [33]. 
With  this  transformation,  the  stress-intensity  factor  can  then  be  calculated 
from  the  following  equation: 

n 

K  =  <7o  y/OZ  y  Ht  F(Ci)  (23) 

i=l 

An  evaluation  was  made  on  the  calculation  of  the  stress-intensity  factor 
using  Equation  (23).  It  was  found  that,  for  n  =  9,  the  computed  K  values 
using  Gaussian  integration  are  essentially  the  same  as  the  ones  obtained  using 
trapezoidal  rule  numerical  integration  with  400  equal  integration  intervals. 


Therefore,  the  more  economical  Gaussian  integration  was  exclusively  used  in 
the  K  calculation  using  the  weight  function  approach. 

An  important  point  on  the  development  of  the  Green's  function  for  attach¬ 
ment  lugs  is  to  be  noted  at  this  point.  In  the  case  of  an  attachment  lug 
loaded  by  a  pin,  which  is  a  contact  problem,  the  pin  bearing  pressure  distri¬ 
bution  between  the  lug  and  the  pin  varies  with  the  crack  length  (Figures 
3-14  through  3-16).  The  pin  bearing  pressure  distribution  has  a  significant 
effect  on  the  unflawed  stress  distribution  (on  the  prospective  crack  surface) 
and  on  the  stress-intensity  factors  (Figure  3-17).  In  the  case  of  finite 
element  analysis  of  cracked  lugs,  the  pin  bearing  pressure  distribution  was 
taken  into  account  in  the  analysis  itself.  In  the  case  of  the  Green's  function 
method,  the  pin  bearing  pressure  distribution  variation  can  be  accounted  for 
by  two  methods. 

The  first  method  is  to  use  the  Green's  function  developed  above  which 
should  strictly  be  used  only  in  conjunction  with  the  correct  unflawed  stress 
distribution  on  the  prospective  crack  surface  for  the  varying  pin  bearing  stress 
distribution  as  the  crack  length  changes.  r<’he  problem  of  calculating  the  pin 
bearing  pressure  distribution  as  the  crack  length  changes  is  statically  in¬ 
determinate,  unless  a  fracture  analysis,  such  as  the  finite  element  solution 
described  before,  is  carried  out.  But,  once  a  fracture  analysis  is  made, 
then  there  is  no  need  for  the  varying  pin  bearing  pressure  distribution, 
because  the  needed  stress-intensity  factors  can  be  obtained  from  the  fracture 
analysis.  Thus,  this  method  becomes  redundant  and  is  discarded. 

The  second  method  is  to  use  available  data  of  stress-intensity  factors 
(from  the  finite  element  procedure,  Figure  3-18)  and  the  unflawed  stress 
distribution  on  the  prospective  crack  surface  (Figure  3-4)  to  modify  the 
Green's  functions  to  account  for  the  varying  pin  bearing  pressure  distribu¬ 
tion.  Such  modifications  have  been  made  lor  the  Green's  function  and  are  pre¬ 
sented  in  Figures  3-26  through  3-28  for  outer-to-inner  radius  ratios  of 
1.50,  2.25  and  3.00,  respectively,  whereas  Figures  3-23  through  3-25  corre¬ 
spond  to  the  original  or  unmodified  Green's  function.  For  calculating  the 
stress-intensity  factors  in  simple  attachment  lugs  loaded  by  oins,  the 
Green's  functions  in  Figures  3-26  through  S-28  should  be  used.  The  Green's 
functions  given  in  Figures  i — 2  3  through  3-25  are  still  useful  and  needed  in 
the  anal  vs  Is  of  attachment  lugs  where  the  unflawed  stress  distribution  on  the 
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prospective  crack  surface  does  not  depend  on  crack  length,  for  example, 
residual  stresses  due  to  the  installation  of  an  interference-fit  bushing. 
Tables  3-4  through  3-9  provide  the  original  and  modified  Green's  function 
values  in  a  tabular  form  at  several  (x-R^/c  locations  and  different  crack 
lengths  for  outer-to-inner  radius  ratios  of  1.50,  2.25  and  3.00. 

In  order  to  put  the  results  in  perspective,  Figure  3-29  compares  the 
stress-intensity  factors  obtained  by  the  compounding  method,  finite  element 
method  and  original  and  modified  Green’s  function  methods.  As  can  be  seen 
from  the  figure,  the  results  obtained  by  the  modified  Green's  function 
method  match  the  finite  element  solutions.  The  results  of  the  original  or 
unmodified  Green's  function  method  tend  to  be  higher,  especially  for  the  R^/R 
ratio  of  1.50.  The  difference  between  the  two  Green's  function  methods 
essentially  reflects  the  effect  of  pin  bearing  pressure  distribution  as  the 
crack  length  changes. 

Figure  3-29  also  compares  the  compounding  method  solution  to  the  other 
solutions.  For  R  /R^  ratios  of  2.25  and  3.00,  the  compounding  solution 
gives  excellent  estimates  of  stress-intensity  factors,  except  when  the  crack 
length  is  very  small.  The  error  at  small  crack  lengths  could  be  corrected 
by  making  proper  use  of  the  stress  concentration  factor  values  for  the  lugs, 
Ktb>  given  in  Figure  3-5.  However,  for  K^R^  =1.50  the  compounding  solu¬ 
tion  tends  to  overestimate  the  stress-intensity  factors,  like  the  original 
Green's  function  method.  This  error  may  be  partially  due  to  the  major  dif¬ 
ferences  between  the  true  pin  bearing  pressure  distribution  and  the  ideal¬ 
ized  distribution  used  in  the  compounding  method. 

2.4  ANALYSIS  OF  LUGS  WITH  INTERFERENCE-FIT  BUSHINGS 

In  order  to  improve  the  crack  initiation  and  crack  growth  life,  the 
concept  of  installing  an  interference- f it  bushing  to  introduce  beneficial 
resicual  stresses  around  the  hole  of  the  lug  prior  to  pin  fitting  has  been 
used  in  aircraft  attachment  lug  design.  For  a  given  fatigue  load  cycle, 
the  installation  of  an  interference  bushing  can  reduce  the  effective  tangen 
tial  stress  range  at  the  likely  location  of  crack  initiation,  resulting  in 
improvement  in  fatigue  and  crack  growth  life.  It  also  reduces  fretting 
damage  of  the  hole  wall  of  the  lug. 
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E  ELEMENT  METHOD 


POUNDING  METHOO 

4FIE0  WEIGHT  FUNCTION  METHOD 

INAL  WEIGHT  FUNCTION  METHOO 


The  stress-intensity  factors  for  a  crack  emanating  from  an  attachment 
lug  having  an  interference-fit  bushing  can  be  estimated  using  the  weight 
function  approach  similar  to  the  one  reported  in  Reference  [34].  The  proce¬ 
dure  consists  of  two  major  steps.  First,  the  effective  unflawed  stress  dis¬ 
tribution  on  the  prospective  crack  surface  was  obtained  by  superimposing  the 
residual  hoop  stresses,  crre>  due  to  the  installation  of  an  intereference-f it 

bushing  on  the  applied  tangential  stresses,  <j  ,  obtained  due  to  the  appli- 

ap 

cation  of  a  pin  loading.  Next,  a  crack  was  introduced  in  this  stress  field 
by  removing  the  tractions  on  the  crack  faces  and  computing  the  corresponding 
effective  stress-intensity  factor  using  the  developed  weight  function,  G, 
and  the  following  equation: 


K  = 


+  arg)  •  G(c,ij)d£ 


(24) 


Unflawed  Stress  Analysis 

The  installation  of  an  interference-fit  bushing  creates  a  compressive 
radial  stress  and  a  tensile  hoop  stress  in  the  lug  in  a  manner  similar  to  a 
thief  wall  cylinder  under  internal  pressure.  On  this  basis,  an  approach 
similar  to  that  of  Seely  and  Smith  [35]  for  a  thick-wall  cylinder  under 
internal  pressure  is  used  to  compute  the  residual  stresses  in  the  lug  due 
to  the  installation  of  an  interference-fit  bushing.  Before  the  residual 
stresses  in  the  lug  can  be  computed,  the  pressure,  P  ,  on  the  surface  of 
contact  between  the  bushing  and  the  lug  must  be  determined. 

Let  R  and  K.  denote  outer  and  inner  radii  of  the  lug  and  r  and 
o  i  °  o 

r^  denote  outer  and  inner  radii  of  the  bushing  before  the  interference-fit 
installation,  respectively.  Let  be  the  inner  radius  of  the  lug  (and 

also  che  outer  radius  of  the  bushing)  after  installation,  and  let  5  be 
the  dltterence  in  these  radii  before  installation. 


5 


K 


r 

o 


(?5) 


After  the  installation  ol  the  bushing,  the  inner  radius  of  the  lug  will 
be  larger  than  its  initial  value  by  an  amount ,  which  is  related  to  the 
unknown  contact  pressure  as 


(26) 


6 


1 


P  R. 
s  1 


where  E  and  jLl  are  the  Young's  modulus  and  Poisson's  ratio  of  the  lug. 

U  L 

At  the  same  time,  the  outer  radius  of  the  bushing  will  be  changed  (decreased) 
by  an  amount,  §2>  which  is  given  by  the  following  equation: 


(27) 


where  E  and  j '1  are  the  Yojng's  modulus  and  Poisson's  ratio  of  the  bush- 

D  D 

ing. 


At  the  completion  of  the  installation  process,  the  difference  in  these 

original  radii,  6dj  disappears  as  a  result  of  the  changes  in  the  length  5. 

K  X 

and  62*  Therefore,  the  sum  of  the  magnitudes  of  5^  and  62  is  equal  to 
5r  ‘-e-  • 


The  contact  pressure  is  obtained  from  the  above  equation  as 


(29) 


With  contact  pressure  determined,  the  residual  stresses  in  the  lug  due 
to  the  installation  of  an  interference-fit  bushing  can  be  calculated  using 
the  following  equations: 


and 


(30) 


P 

s  1 

2  2 

R  -  R, 
o  1 


P  K1 
s  1 

2  2 

R  -  R, 
o  1 


x  >  R^ 


where  x  Is  the  radial  distance  as  shown  in  Figure  3-30. 

Equation  (30)  can  be  used  to  compute  the  residual  tangential 
stresses  along  the  radial  direction  (x-axis)  in  attachment  lugs  due  to  the 
installation  of  interference-fit  bushings.  In  the  following  discussion, 
results  are  presented  for  attachment  lugs  with  lug  outer  radius  to  bushing 
inner  radius  ratio  (R^r^)  of  2.25.  The  nondimensional  parameters  con¬ 
sidered  in  the  analysis  are  interference  levels  (6u/r.)>  bushing-to-lug 
rigidity  (E^/E,  )  and  bushing  thickness  t^/r^).  The  effects  of  these  param¬ 
eters  on  the  stress  distribution  along  the  x-axis  in  attachment  lugs  are 
presented  In  Figures  3-30  through  3-32.  Figure  3-30  shows  the  Increase  in 
residual  stresses  with  increasing  diametral  interference  for  a  bushing-to- 
lug  rigidity  ratio  (E  /E,  )  of  3.0  and  a  bushing  thickness  (t_/r  )  of  0.12. 
Figure  3-31  shows  the  increase  in  residual  stresses  with  incrasing  bushing- 
to-lug  rigidity  ratio  ( E  / E  )  for  a  bushing  thickness  (t  /r  )  of  0.12  and 

an  interference  level  (6.,/r.)  of  0.00533.  Figure  3-32  shows  the  increase 

K  1 

in  residual  stresses  with  increasing  bushing  thickness  (t  /r.)  for  a 

d  i 

bushing-to-lug  rigidity  ratio  (E  /E  )  of  3.0  and  an  interference  level 

D  l* 

(5  Jr  )  of  0.00533. 

K  i 

The  tangential  stresses  along  the  x-axis  due  to  the  application  of  pin 
loading  are  then  computed  using  the  finite  element  method.  In  the  analysis, 
it  is  assumed  that  the  bushing  and  the  lug  remain  in  contact  and  that  no 
slippage  occurs  along  the  hole  wall  surface  during  the  application  of  the 
load.  The  computed  tangential  stress  normalized  by  the  far-field  applied 
stress  for  a  lug  with  a  of  2.25  are  shown  in  Figures  3-33  and  3-34. 

Figure  3-11  depicts  the  decrease  In  stresses  due  to  pin  loading  with  In¬ 
creasing  bushing-to-lug  rigidity  ratio  (Eu/E, )  tor  a  constant  bushing  thlek- 
ness  (t^/r^)  af  0.17.  Figure  3-34  presents  the  decrease  In  stresses  due  to 

pin  loading  with  increasing  bushing  thickness  (t  /r.)  for  a  constant  bushing- 

11  i 

to-lug  rigidity  ratio  (Eu/E, )  ot  3.0 

b  i. 
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NON-DIMENSIONAL  RESIDUAL  TANGENTIAL  STRESS,  <*y/E 
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Figure  l-U.  Stress  Distributions  Along  x-Axis  in  an  Attachment  tug 
With  lushing  tor  Various  Su»hing-lo-Lug  Rigidity  Ratios 


(x-R,)/(R0-R,) 


Figure  3-33.  Stress  Distributions  Along  x-Axis  in  an  Attachment  Lug 
With  Bushing  Due  to  Pin  Loading  for  Various 
Bushing-to-Lug  Rigidity  Ratios 


The  results  presented  in  Figures  3-30  through  3-34  pertaining  to  the 
stress  analysis  of  lugs  with  bushings  should  be  used  with  caution  in  the 
sense  that  they  can  be  used  as  long  as  the  assumptions  are  not  violated. 

In  other  words,  these  results  can  be  used  as  long  as  there  is  sufficient 
interference  to  prevent  separation  between  the  bushing  and  the  lug.  When 
the  interference  level  is  small  and/or  when  the  crack  length  is  large  there 
will  be  separation,  which  again  depends  on  the  geometry  and  the  magnitude  of 
the  applied  loading.  Also,  the  above  results  cannot  be  used  when  the  total 
effective  stress  (residual  stress  plus  applied  stress)  exceeds  the 
material  yield  strength.  If  the  total  effective  stress  exceeds  the  material 
yield  strength,  a  nonlinear  elasto-plastic  finite  element  analysis  such  as 
used  in  Reference  [34]  should  be  conducted  to  obtain  the  total  effective 
stress.  In  the  current  study,  only  the  linear  analysis  is  performed.  Once 
the  total  tangential  stresses  are  obtained  by  superposing  residual  stresses 
and  the  applied  stresses,  then  they  can  be  used  to  calculate  the  stress  in¬ 
tensity  factors  for  a  crack  emanating  from  an  attachment  lug  having  an 
interference-fit  bushing  by  using  Equation  (24)  where  G(c,£)  is  the  origi¬ 
nal  or  unmodified  Green's  function. 

Stress-Intensity  Factor  Analysis 

To  study  the  effects  of  various  parameters  on  actual  (dimensional) 

stress-intensity  factors  rather  than  nondimensional ,  a  physical  lug 

having  an  outer  radius  of  1.6875  inches  with  a  bushing  inner  radius  of 

0.75  inch  (R  /r,  =  2.25)  subjected  to  a  far  field  stress  (a)  of  6  ksi  and 
o  i  o 

stress  ratio  (Kc  )  of  0.1  is  considered  in  this  discussion.  For  this  lug 
tar 

configuration  and  loading,  an  effective  stress  ratio  (K^^)  at  the  lug  hole 

wall  (x  =  Rj)  as  a  function  of  interference  levels  can  be  computed  and  is 

presented  in  Figure  3-35a.  Figure  3-35b  shows  the  variation  of  Reff  along 

the  x-axis  of  the  lug  for  an  interference  level  (5^/r^)  of  0.00533.  It  is 

seen  from  the  figure  that  R  , ,  is  constant  almost  throughout  the  net  sec— 

eft 

tion  (x-axis)  except  near  the  outer  surface  of  the  lug. 

The  stress-intensity  factors  were  computed  using  Equation  (24)  and  are 

presented  in  Figures  3-36  and  3-37.  Figure  3-36  shows  the  effective  stress- 

intensity  factor  range,  AK,  for  various  bushing-to-lug  rigidities  (Eg/E^) 

for  a  given  bushing  thickness  (t^/r^)  of  C.12.  The  data  corresponding  to 

ED/E,  =  0  represent  a  simple  lug  with  no  bushing.  It  is  clear  from 
B  U 
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to 


Figure  3-35.  Effective  Stress  Katio  in  an  Attachment  Lug  with  Bushing 
for  Various  Interference  Levels  and  Along  x-Axis 
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STRESS  INTENSITY  FACTOR  RANGE  AK,  KSI  VlN 


5 


Figures  3-35a  and  3-36  that  the  installation  of  the  interference-fit  bushing 
causes  an  increase  in  the  effective  stress-intensity  factor  ratio,  = 

<3  ,  to  =  K  ,  /K  ,  but  a  significant  decrease  in  the  stress- intensity 
factor  range, AK.  This  will  result  in  the  reduction  of  fatigue  crack  growth 
rate.  The  effects  of  bushing  thickness  (t_/r.)  on  the  effective  stress- 
intensity  factor  range,  AK,  are  presented  in  Figure  3-37.  Based  on  the  re¬ 
sults  shown  in  Figures  3-36  and  3-37,  it  may  be  concluded  that:  (1)  for 
constant  bushing  thickness,  an  increase  in  the  bushing  rigidity  decreases  the 
effective  stress-intensity  factor  range;  and  (2)  for  constant  bushing  rigid¬ 
ity,  an  increase  in  bushing  thickness  decreases  the  effective  stress-intensity 
factor  range. 

3.  STRESS-INTENSITY  FACTORS  FOR  CORNER  CRACKS 

One  of  the  most  common  types  of  flaws  for  which  there  exists  no  closed 
form  analytical  solution  is  the  corner  crack  at  a  circular  hole.  To  date, 
several  approximate  methods  have  been  proposed  for  computing  the  stress- 
intensity  factors  of  a  quarter-elliptical  crack  emanating  from  the  comer  of 
an  open  hole  located  in  the  geometric  center  of  a  plate.  The  methods  range 
from  an  empirical  equation  which  was  developed  using  the  fatigue  crack 
growth  method  of  calibrating  the  measured  crack  growth  rate,  da/dN,  and  the 
stress-intensity  factor  range,  AK,  to  one-  and  two-dimensional  compounded 
solutions,  to  tho  sophisticated  three-dimensional  finite  element  analysis. 

In  this  section,  three  methods  in  an  increasing  level  of  complexity  and 
sophistication  are  presented  for  the  determination  of  the  stress-intensity 
factors  for  single  corner  cracks  in  aircraft  attachment  lugs. 

3 A  ONK-PARAMKTKR  COMPOUNDING  APPROXIMATION 

The  corner  crack  is  more  difficult  to  analyze  than  the  through-the- 
thickness  crack  because  it  has  both  surface  length  and  bore  depth  dimensions 
and  both  must  be  considered  in  the  crack  growth  analysis.  0,.\  possible 
simplification  is  to  assume  a  fixed  relationship  between  crack  depth  "a"  and 
surface  length  "c",  so  that  only  one  crack  length  parameter  is  independent. 
This  was  the  approach  used  in  deference  [ 1 4 ]  to  analyze  corner  cracks  at 
fastener  holes.  The  following  stress-intensity  correction  factor,  governing 
the  growth  of  the  corner  crack  in  the  length  direction,  was  verified  empir¬ 
ically  on  open  hole  fatigue  coupons  in  Reference  [14]: 


0.2886 


(31) 


(c) 

<i>  =  i- 

71  1  *  2<f>2<|)2 

This  factor  can  be  multiplied  by  the  through-the-thickness  crack  solu¬ 
tion,  °btafn«d  by  the  compounding  method  and  given  in  Equation  (10). 

k(c>  "  \uc  *n  <”> 

Equation  (32)  provides  a  simple  one-parameter  K  estimate  for  the  point 
near  the  lug  face.  This  factor  can  be  used  even  after  a  >B,  to  smoothly 
change  K  as  the  comer  crack  becomes  a  through-the-thickness  crack. 

Equation  (32)  was  used  to  estimate  the  stress-intensity  factors  for 

corner  cracks  at  attachment  lugs  having  R  /R,  ratios  of  1.50,  2.25,  and 

o  l 

3.00.  These  results  are  tabulated  In  Tables  3-10  and  3-11  for  B/ft ^  =  2/3 
and  i/3,  respectively.  Note  that,  in  the  calculation,  the  depth-to-length 
ratio  (a/c)  was  assumed  to  be  constant  and  equal  to  1.33,  which  was  found 
exDerimentally  by  Schijve  and  Hoeymakers  [2]. 

3.2  TWO-PARAMETER  WEIGHT' FUNCTION  APPROXIMATION 

In  an  earlier  contractual  program  with  the  Flight  Dynamics  Laboratory,  Hsu 
et  a l  [36]  developed  a  simple  procedure  for  estimating  the  stress-intensity 
factors  along  the  boundary  of  quarter-elliptical  corner  cracks  emanating 
from  fastener  holes  from  corresponding  through  crack  solutions.  Correlation 
between  the  calculated  stress-intensity  factors  and  those  deduced  from  the 
tests  were  good.  A  similar  procedure  can  be  used  to  estimate  the  stress- 
intensity  factors  along  the  periphery  of  a  quarter-elliptical  corner  crack 
in  an  attachment  lug. 

The  stress-intensity  factor  for  a  single  through-the-thickness  crack 
of  length  c  in  an  attachment  lug  is  given  by 

K(jr~)  =  crQ  >/5tc  (jp) 

where  (c/kj)  is  the  through-the-thickness  correction  factor. 

Then  the  corresponding  stress-intensity  factor  for  a  single  quarter- 

elliptical  corner  crack  can  be  written  as 
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for  —  <  1 
c 


and 


9t(^r)  x  (Hm4v 


(35) 


for  -  >  1 
c 


where 


x-R 

M  -  F  +  - -  (1-F) 

c  c 


(36) 


in  which  a  is  the  elliptical  angle  measured  from  the  hole  wall,  x  is  the 

distance  from  the  center  of  the  hole  to  the  particular  point  of  interest  on 

the  crack  periphery,  M|/0  is  the  combined  front  free  surface  and  flaw  shape 

factor,  M  is  the  wall  curvature  correction  factor,  and  is  the  back 
c  D 

surface  correction  factor.  The  values  of  Mj  /<P  andF,  and  are  presented  in 
References  [36]  and  [58]  ,  respectively.  As  in  Reference  [36]  ,  it  is 
assumed  that  for  a  given  number  of  applied  load  cycles,  the  extension  of  the 
quarter-elliptical  crack  border  is  controlled  by  the  stress  intensity  factors 
at  two  points;  namelv.  the  intersections  of  the  crack  periphery  with  both  the 
hole  wall  and  plate  surface  (i.e.,  and  K^.).  In  general,  the  stress 
intensity  factors  at  these  two  locations  are  different,  resulting  in  differ¬ 
ent  crack  growth  rates.  Therefore,  the  new  flaw  shape  aspect  ratio  after 
each  crack  growth  increment  differs  from  the  preceding  one.  The  new  flaw 
shapt  aspect  ratio  is  computed  using  the  new  crack  lengths  at  both  the  hole 
wall  and  plate  surface.  All  the  correction  factors  have  been  reduced  to 
simple  equations  for  locations  A  and  C,  anti  the  two-parameter  corner  crack 
stress  Intensity  factot  formulas  arc  summarized  in  Figure  3-38. 


in  Reference  fit*!,  the  above  process  was  repeated  until  the  crack  length 
along  the  hole  wall  was  equal  to  the  plate  thickness.  At  that  time,  the 
craek  was  assumed  to  be  a  through-the-thickness  crack  with  length  c.  This 
assumption  was  based  upon  the  experimental  observation  that  after  the 
craek  penetrates  the  back  surface  and  the  cyclic  load  application  continues, 
the  back  surface  crack  length  Increases  much  faster  than  that  of  the  front 
sutface  until  the  front  of  the  through-the-thickness  crack  becomes  stable. 
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Figure  3-38.  Two-Farameter  Corner  CracV  Stress  Intensity  Factor  Formulas 


This  transitional  crack  growth  assumption  is  very  reasonable  for  a  crack 
at  a  fastener  hole,  since,  at  the  time  the  crack  grows  through  the  thickness, 
the  remaining  net  section  is  usually  much  larger  than  the  surface  crack 
length.  However,  for  an  atta-  lunent  lug.  when  a  part -through  crack  grows 
through  the  thickness,  the  remaining  net  section  is  usually  small  in 
comparison  to  the  surface  crack  size.  Therefore,  a  proper  transitional  crack 
growth  criterion  Is  needed  for  the  transition  pev iod  from  the  time  the  crack 
penetrates  the  back  surface  to  the  time  the  crack  lengths  are  essentially 
equal  on  the  front  and  back  surfaces. 


Figure  3-39  shows  the  transitional  crock  geometry,  in  which  c^  and 
c^  are  crack  lengths  on  the  front  and  back  surfaces,  respectively. 
Collipriest  and  Ehret  [37]  proposed  a  stress-intensity  magnification  factor 
for  the  crack  tip  at  the  back  surface  of  a  surface  crack  as 


a 

t 


i 


1/2 


2 


for  c_  >  0 
B 


(37) 


When  the  back-side  crack  length  equals  the  front-side  crack  length,  the 
magnification  is  unity  and  the  through  crack  has  achieved  a  uniform  front. 
This  factor  will  be  used  to  estimate  the  stress-intensity  factor  at  the  back 
surface  of  the  transitional  crack,  i.e. 


“o  -  »iT  <^>  •  <1, 


(38) 


The  st ress-intensity  factor  at  the  front  surface  will  be  calculated 
using  Equation  ( 19)  and  an  imaginary  crack  length  along  the  hole  wall  direc¬ 
tion,  a  * ,  i.e. 


\  3"er 


c  H*  c 

(i-  (~>  ~  <~.  W°) 

T  4>  a 


(3*)) 


An  imaginary  crack  length  a'  can  be  determined  by  fitting  an  ellip¬ 
tical  equation  through  points  C  and  B,  as 
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-1/2 


a’ 


-r  n 


(40) 


When  the  back-side  crack  length,  cn,  approaches  that  of  the  front  side  (i.e. 

D 

B"C"  as  shown  in  Figure  3-38),  the  factor  c^/a'  becomes  very  small  and  the 
corresponding  M  j  /tf>  ratio  approaches  unity.  After  that,  the  transitional 
crack  becomes  a  through-the-thickness  crack  with  a  uniform  front,  i.e. 

<41) 

Equations  (34)  and  (35)  were  used  to  compute  the  stress-intensity 
factors  at  tne  intersections  of  the  crack  periphery  and  the  hole  wall  and 
plate  surfaces.  The  2-D  finite  element  method  solution  results  from  Table  3-3 
were  used  for  the  basic  through-the-thickness  crack  solution,  K/c/R^).  The 
results  are  shown  in  Figures  3-40  through  3-42  for  R  /R^  ratios  of  1.50, 

2.25  and  3.0,  respectively.  In  each  plot,  the  flaw  shape  aspect  ratio  a/c 
was  assumed  as  constant.  It  should  be  noted  that  the  computed  part-through 
ct  >.ck  stress-intensity  factors  are  normalized  by  cf^  *  These  solutions 

correspond  to  corner  crack  solutions  for  a<B.  Figure  3-43  presents  the 
normalized  stress-intensity  factors  at  the  lug  surface  having  a  constant 
a/c  ratio  of  1.33  and  a  KQ/l<i  ratio  ranging  from  1.5  to  3.0.  In  this 
figure,  the  solutions  are  presented  even  after  the  crack  breaks  through  the 
back  side,  i.e.  a  > B.  Corresponding  normalized  stress-intensity  '  factors 
computed  using  the  one-parameter  coinpounding  method  are  also  shown  in  the 
figure  for  comparison. 


3.3  't'HRKF.-tHMENS  I0NA1.  CRACKED  FINITE  ELEMENT  METHOD 

Numerical  methods  such  as  the  slicing  technique  [38],  the  conventional 
three-dimensional  (3-U)  finite  element  method  [39],  the  boundary  integral 
equation  approach  [4ol,  and  the  3-D  alternating  technique  [41 ]  have  been 
used  to  analyze  corner-crack  problems.  To  accurately  depict  the  extreme 
stress  gradient  existing  in  the  vicinity  of  a  crack  tip,  the  necessary  com¬ 
puter  effort  in  each  case  is  considerable.  In  an  attempt  to  minimize  the 
economic  problem  of  the  conventional  approach  and  to  extract  the  stress- 
intensity  (actors  directly,  a  3-dimensional  cracked  finite  element  was  used 
which  is  capable  of  characterizing  the  crack-tip  stress  singularity  intern¬ 
ally.  The  procedure  used  in  the  development  of  the  3  D  cracked  element  is 
described  m  keierence  [42]  and  summarized  as  follows. 
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Figure  3-40.  Normalized  Stress  Intensity  Factors  for  Corner  Cracks 

in  a  Straight  Attachment  Lug  Having  a  R  /R.  Ratio  of  1.50 
and  a  B/R^  Ratio  of  2/3  ° 


Two  types  of  20-node  3-D  cracked  elements,  as  shown  in  Figure  3-44, 
were  developed  using  i  hybrid  displacement  finite  element  procedure.  The 
variational  principle  which  governs  the  assumed  displacement  hybrid  finite 
element  model  is  the  stationary  condition  of  a  modified  total  potential 
energy  functional  with  a  relaxed  requirement  of  interelement  boundary  dis¬ 
placement  continuity,  a  priori.  This  variational  principle  is  a  three-field 
variable  principle.  The  three  field  variables  are  the  element  interior  dis¬ 
placements,  element  boundary  displacements,  and  the  Lagrange  multipliers. 

The  Lagrange  multipliers  are  physically  the  interelement  boundary  tractions 
and  are  assumed  to  match  the  independently  assumed  element  interior  and 
boundary  displacements  at  the  interelement  boundary.  Three-dimensional 
asymptotic  solutions  for  displacements  (\/ r  -  type)  and  stresses  (l / \/  r  - 
type)  near  the  crack  front  are  embedded  in  this  procedure.  Apart  from  satis¬ 
fying  relevant  field  equations,  the  variational  principle  also  enforces  the 
conditions  of  displacement  continuity  and  traction  reciprocity  at  the  inter¬ 
element  boundary  a  posteriori,  assuring  the  convergence  of  the  finite  ele¬ 
ment  procedure.  The  procedure  is  capable  of  analyzing  mixed  mode  fracture 
problems  (Modes  l,  II  and  III),  and  the  three  stress-intensity  factors  at 
various  locations  along  the  crack  front  are  treated  as  unknowns  along  with 
the  generalized  nodai  displacements  of  the  structure.  The  final  set  of 
algebraic  equations  governing  the  global  nodal  displacements  and  the  three 
stress-intensity  factors  at  various  locations  along  the  crack  front  can  be 
written  as  follows: 

[Kj]  M  *  [k/  |K*|  »  {Qjj  (42) 

[K2]  |q*|  *  ]  | K* j  =  |Q?|  (43) 

Where  |q*j  are  the  structure's  global  nodal  displacements,  {K*|  are 
the  mixed  mode  stress-inter. !  ■  '»rcrs  at  various  locations  along  the  crack 
front,  [K  ],  [K,]  ami  [fC^]  are  the  corresponding  stiffness  matrices  (super¬ 
script  T  represents  the  transpose),  and  JqJ  and  {Q?}  are  the  corresponding 
nodal  forees.  It  is  evident  from  Equations  (42)  and  (43)  that  the  solution 
for  stress-intensity  factors  can  be  obtained  directly  from  the  finite  ele¬ 
ment  solution  procedure.  This  eliminates  additional  post-processing  of  dis¬ 
placement  or  stress  solutions  to  obtain  the  stress-intensity  factors  through 
methods  such  as  the  crack  opening  displacement  method,  nodal  force  method, 
ete.  a  detailed  description  of  the  hybrid  displacement  procedure  was  docu¬ 
ment  ed  in  keference  [42], 
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Figure  3-44.  20-Node  Yhtee-Diwsns tonal  Cracked  IClewmts 


The  accuracy  and  convergence  of  the  three-dimensional  hybrid  displace¬ 
ment  finite  element  procedure  have  been  tested  and  verified  through  the  solu¬ 
tion  of  several  complex  fracture  problems  of  interest  in  aerospace  structural 
components,  nuclear  pressure  vessels  and  components,  and  solid  rocket  motor 
grain  applications  [43-45].  Raju  and  Newman  [46],  using  Tracey's  [39]  wedge 
shaped  distorted  isoparametric  crack  elements,  made  a  convergence  study  for 
the  solution  of  semielliptical  surface  flaws  in  thin  plates.  A  similar  con¬ 
vergence  study  was  also  made  using  the  present  three-dimensional  hybrid 
finite  element  procedure  in  Reference  [44].  As  discussed  in  Reference  [44],  the 
present  solution  with  4555  degrees  of  freedom  has  excellent  agreement  with  that 
of  Raju  and  Newman  for  their  highest  degrees  of  freedom,  6867. 

The  three-dimensional  cracked  element  was  used  to  compute  the  stress- 
intensity  factors  for  single  quarter-elliptical  corner  cracks  at  straight 
lugs.  A  total  of  eight  problems  of  corner  cracks  in  attachment  lugs  with 
various  outer-to-inner  radius  ratios,  crack  aspect  ratios,  and  crack  depth- 
to-lug  thickness  ratios  were  considered  for  the  present  analysis.  The  geom¬ 
etry  of  an  attachment  lug  with  a  corner  crack  is  given  in  Figure  3-45 (a),  and 
the  various  parameters  for  the  eight  problems  considered  are  defined  in 
Table  3-12.  The  bearing  pressure  at  the  lug  hole  wall  was  assumed  to  be 
uniform  and  acts  only  on  the  right  half  of  the  hole  as  shown  in  Figure  3-45(a). 
This  assumption  of  uniform  bearing  pressure  over  the  right  half  of  the  lug 
hole  wall  was  made  due  to  the  complexity  involved  in  the  generation  of  equiv¬ 
alent  nodal  forces  c  esponding  to  an  actual  pin-lug  contact  bearing  pres¬ 
sure  distribution  a  three-dimensional  case.  Consequently,  the  solutions 
by  the  two-dimensional  Green's  function  method  generated  for  comparison 
purposes  also  correspond  to  the  same  bearing  pressure.  Also,  the  original 
Green's  functions  were  used  in  the  computation  of  two-dimensional  solutions. 

3'he  assumption  of  same  bearing  pressure  in  both  cases  is  important,  because 
the  stress-intensity  factor  solutions  vary  significantly  depending  upon  the 
bearing  pressure  distribution,  namely,  uniform,  cosine  and  actual  pin-lug 
pressure  d ist r ibut ions . 

A  typical  finite  element  breakdown,  number  of  elements  and  total  number 
of  degrees  of  freedom  for  the  present  analysis  are  presented  in  Figure  3-45  (b)  . 

In  Figures  3-46  through  3-50,  solutions  of  stress-intensity  factor  variations 
along  the  crack  front  by  both  two-  and  three-dimensional  procedures  are  pre¬ 
sented  tor  problem  numbers  1  and  2,  3,  4  and  5,  6  and  7,  and  8,  respectively. 
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(b)  Finite  Element  Breakdown 


.  Geometry  and  Finite  Element  Breakdown  ot 
with  a  Corner  Crack  for  Three-Dimensional 
Element  Analysis 


an  Attachment 
Cracked  Finite 
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TABLE  3-12.  CORNER  CRACK  PROBLEMS  ANALYZED 
BY  THREE-DIMENSIONAL  CRACKED 
FINITE  ELEMENT  PROCEDURE 


PROBLEM 

NUMBER 

1 

2 

3 

4 

5 

6 

7 

8 


Vl 

a/c 

'  . . ■  — 

a/B 

1.50 

1.0 

0.50 

1.50 

1.5 

0.50 

2.25 

1.0 

0.1667 

2.25 

1.0 

0.25 

2.25 

1.5 

0.25 

2.25 

1.0 

0.50 

2.25 

1.5 

0.50 

2.25 

1.5 

0.75 

N. 

N, 


~  -  0.1667 

D 


ELLIPTICAL  ANGLE,  a 


lt,ure  1-4'/.  Nor-nallzcd  Seres*  Intensity  Factors  for  Corner  Cracks 

at  a  Straight  Attachment  hug  Having  a  R  /R  Ratio  of  2. 25 
and  an  */B  Ratio  of  0.1667  ° 


ELLIPTICAL  ANGLE,  Cl 


Normalized  Stress  Intensity  Factors  for  Corner  Cracks 
at  a  Straight  Attachment  Lug  Having  a  K  /R.  Ratio  of  2.25 
and  an  a/B  Ratio  of  0.25  ° 


The  stress-intensity  factor  solutions  are  normalized  by  the  exact  solution 
of  a  through-the-thickness  crack  of  length  c  in  an  infinite  plate  subjected 
to  a  crack  pressure  of  <7  .  Elliptical  angles  (  0 i  )  of  0°  and  90°  refer 
to  the  intersections  of  the  crack  front  with  the  hole  wall  and  lug  surface, 
respectively . 

A  comparison  of  the  results  in  Figures  3-46  through  3-50  reveals  that  the 
two  solution  procedures  agree  well  at  d  =  0°,  but  significant  differences 
exist  near  Ct  =  90°.  For  lugs  with  an  outer-to-inner  radius  ratio  of  1.5, 
the  maximum  difference  at  Ot  =  90°  is  about  42  percent  and  the  difference 
decreases  as  the  aspect  ratio  of  the  crack  (a/c)  increases.  For  lugs  with 
outer-to-inner  radius  ratio  of  2.25,  the  maximum  difference  at  01=90°  is 
about  26  percent  and  the  difference  decreases  as  the  aspect  ratio  of  the 
crack  (a/c)  and  the  crack  depth-to-thickness  ratio  (a/B)  increase.  The  overall 
behavior  of  the  differences  in  the  solutions  is  that  they  decrease  as  the 
outer-to-inner  radius  ratio,  crack  aspect  ratio  and  crack  depth-to-thickness 
ratio  increase.  Also,  the  angle  at  which  the  normalized  stress-intensity 
factor  is  minimum  is  lower  for  the  two-dimensional  procedure  in  almost  all 
of  the  problems.  Thus,  from  these  results,  the  two-dimensional  procedure 
seems  to  he  consistently  overestimating  the  stress-intensity  factor  solutions 
except  for  elliptical  angles  close  to  0°. 

The  one-parameter  compounding,  two-parameter  weight  function  and  the 
three-dimensional  cracked  finite  element  solutions  are  compared  in  terms  of 
corner  crack  correction  factors  at  the  lug  surface  in  Figure  3-51.  The  corner 
crack  correction  factors  are  obtained  by  normalizing  Kc  with  the  stress  inten¬ 
sity  factor  (by  weigh,  function  method)  of  a  through-the-thickness  crack  of 
length  c  and  are  expressed  as  a  function  of  a/B. 

4.  ELASTOPLASTIC  ANALYSIS 

Several  analytical  methods  of  varying  complexities  have  been  developed 
and  discussed  in  previous  sections  for  through-the-thickness  and  corner 
cracks  m  straight  attachment  lugs.  These  analyses  are  basically  for 
linear  fracture  mechanics  problems.  However,  attachment  lugs  may 
undergo  substantial  plastic  yielding  around  the  lug  hole,  depending  on  the 
specific  application  and/or  design.  For  example,  the  concepts  of  cold¬ 
working  of  the  hole  or  installation  of  an  interference-fit  bushing  prior  to 
pin  fitting  have  been  extensively  used  in  actual  aircraft  attachment  lug 
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design  practices  to  improve  the  crack  initiation  life.  These  concepts  are 
basically  used  to  introduce  fatigue-improving  residual  stresses  around  the 
hole  of  the  lug.  The  lugs  may  also  undergo  plastic  yielding  due  to  the 
application  of  high  pin  loads.  It  is  important  to  develop  analytical  pro¬ 
cedures  for  analyzing  such  problems. 

The  analytical  procedure  is  similar  to  that  developed  and  verified  for 
through-the-thickness  cracks  emanating  from  interference-fit  and  cold-worked 
fastener  holes  in  infinite  plates,  under  a  previous  Air  Force  contract 
(Reference  [36]).  The  method  involves  an  initial  elasto-plastic  stress  analysis 
of  attachment  lugs  to  obtain  the  stress  distribution  on  the  prospective  crack 
surface.  These  stresses  are  then  used  in  conjunction  with  the  Green's  function 
to  calculate  the  stress  intensity  factors  for  cracks  emanating  from  lug  holes. 
The  Green's  function  developed  for  linear  problems  is  also  used  for  the 
elasto-plastic  case. 

The  elasto-plastic  stress  analysis  has  been  carried  out  for  two  mate¬ 
rials,  4340  steel  and  7075-T651  aluminum,  which  are  being  used  for  testing 
in  the  present  contract  work.  These  materials  have  been  characterized  and 
the  nonlinear  elasto-plastic  stress-strain  relationships  were  approximated 
by  trilinear  representation  and  are  given  in  Figure  3-52.  The  unloading 
moduli  have  been  assumed  to  be  the  same  as  the  initial  elastic  moduli  for 
both  of  the  materials.  The  analysis  was  performed  for  lugs  with  outer-to- 
inner  radius  ratios,  Ro/K^,  °f  1.50,  2.25,  and  3.00  and  a  constant  lug  inner 

radius  of  0.75  inch.  Far-field  maximum  stresses,  (j  ,  of  35  and  15  ksi  were 

o 

selected  for  steel  and  aluminum  lugs,  respectively.  Isotropic  hardening 
with  the  von  Mises  yield  criterion  was  used  in  the  elasto-plastic  analysis. 
Steel  pins  were  used  in  ail  cases,  and  it  was  also  assumed  that  the  pin  did 
not  yield. 

The  finite-element  method  was  again  used  for  the  elasto-plastic 
stress  analysis  of  attachment  lugs.  The  finite  element  breakdown  and 
the  use  of  spring  elements  with  high  stiffness  to  simulate  the  pin- 
lug  contact  surface  were  similar  to  those  in  the  iinear  stress  analysis. 

The  model  was  loaded  by  a  single  concentrated  force  at  the  center  of  the  pin. 

In  the  nonlinear  analysis,  a  loading  was  first  applied  to  reach  the  yield 
point,  followed  by  incremental  loadings  to  reach  the  maximum  tar-field 
stresses  listed  above.  The  loadings  were  then  reduced  to  give  stress  levels 
corresponding  to  far-field  stress  ratios,  R,  of  0.5  and  O.l. 
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(a)  STRESS-STRAIN  CURVE  FOR  4340  STEEL 
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(b)  STRESS-STRAIN  CURVE  FOR  7075-T651  ALUMINUM 


Figure  }-52.  Material  Properties  Used  in 

Elasto-Plastic  Stress  Analysis 
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In  the  plastic  region,  the  following  two  convergence  criteria  were  used 
in  the  analysis.  A  tolerance  criterion  of  AE  /AE  <  0.1  was  used  in 

It  we  1 

the  analysis,  where  AE^,  AE^  ^  are  the  energy  changes  for  the  Iteration 
and  the  load  increment,  respectively.  Regardless  of  the  above  tolerance 
criterion,  the  analysis  was  iterated  at  least  twice  for  each  load  incre¬ 
ment.  In  all  the  solutions  presented  here,  the  maximum  value  of  A  ^ / 

AE  was  less  than  0.035. 

!«•  1 

The  unflawed  stress  distribution  along  the  prospective  crack  surface, 

90°  away  from  the  load  line,  obtained  from  the  elasto-plastlc  finite  element 

analysis  for  the  far  field  stress  levels  given  before  are  presented  in 

Figures  3-53  through  3-55  for  steel  lugs  and  Figures  3-56  through  3-58  for 

aluminum  lugs.  These  results  are  given  in  increasing  order  of  R  /R  ratios 

o  i 

of  1.50,  2.25  and  3.00. 

For  a  steel  lug  with  R  /R,  of  1.50,  the  stresses  (  a  )  corresponding 

o  i  max 

to  maximum  far-field  stress  of  35  ksi  and  the  stresses  (a  ,  R  «  0.5  and 

min 

<Tmin»  R  =  0.1)  corresponding  to  minimum  far-field  stress  for  stress  ratios 

of  0.5  and  0.1,  respectively,  are  given  in  Figure  3-53.  From  these  stresses, 

the  difference  between  maximum  and  minimum  stresses  for  stress  ratios  of  0.5 

and  0.1,  4<7R=Q  5  and  4<fR_0  ^ ,  respectively,  can  be  obtained  and  are  also 

included  in  the  figure.  At  the  maximum  load,  about  20  percent  of  the  lug 

ligament  has  yielded.  The  effects  of  plastic  yielding  of  the  lug  can  also 

be  seen  in  the  distributions  of  <J  ,  _  .  and  O  ,  „  _  While 

min,  R=0.5  min,  R=0.1 

unloading,  due  to  plastic  yielding,  the  stresses  near  the  lug  hole  decrease 

at  a  higher  rate  than  the  rest  of  the  ligament  and  become  negative  at  the 

minimum  stress  corresponding  to  R  *  0.1.  If  the  loading  was  completely 

removed,  significant  negative  residual  stresses  would  exist  near  the  lug 

hole.  In  calculations  of  crack  growth  any  compressive  minimum  stresses  near 

the  lug  hole  are  neglected  due  to  eraek  closure,  and  A<J  ,  is  assumed 

K*0  •  1 

not  to  exeeed  o  .  However,  for  illustration  purposes,  the  aetual 
difference  in  stresses  near  the-  lug  hole  is  indicated  in  Figure  3-53  by  a 
dotted  line.  Corresponding  solutions  for  an  aluminum  lug  with  the  same 
geometry  arc  presented  in  Figure  3-56,  and  the  results  are  qualitatively 
the  same. 

Similar  analytical  solutions  for  lugs  with  an  R  /R  ratio  of  2.25  are 

o  l 

given  in  Figures  3-54  and  3-57  for  steel  and  aluminum  materials,  respec¬ 
tively.  In  this  case,  tor  both  the  materials,  the  lugs  yield  only  slightly: 
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Figure  3-53.  Stress  Distribution  Along  x-Axls 
for  a  Steel  Lug  with  l^/K^  =  1*5° 
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S7.  Stress  Distribution  Along  x-Axis 

for  an  Alualnun  Lug  with  it  /it  a 

o  i 
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that  is,  about  5  percent  of  the  lug  ligament.  Also,  for  R=0.1,  the  minimum 

stress  near  the  lug  hole  does  not  become  negative  for  this  lug  configuration. 

Thus,  Actr_q  ^  and  AaR0  presented  are  the  actual  differences  in  the 

stresses.  Results  for  steel  and  aluminum  lugs  with  a  R  /R,  value  of  3.00  are 

0  1 

presented  in  Figures  3-55  and  3-58,  respectively.  In  this  case,  the  lugs 
barely  yield,  and  very  small  plastic  strains  were  calculated  by  the  elasto- 
plastic  analysis. 

The  above  stress  ranges,  computed  from  the  elasto-plastic  analysis,  are 
then  integrated  with  the  Green’s  function  to  compute  the  stress  intensity 
factors  for  through-the-thickness  cracks  in  lugs  that  are  loaded  above  vield. 

The  stress  intensity  factors  can  be  calculated  using  the  Green's  func¬ 
tion  method  by  the  following  equation: 

K(c)  =  aQ  V^c  J  G(«,$)  ff  ($)  d£  (44) 

a 

The  stress  intensity  factor  range  for  through-the-thickness  cracks  can 
be  computed  from  the  above  equation  by  either  inputting  the  maximum  and 
minimum  stresses  and  subtracting  the  values,  or  more  directly,  by  using  the 
stress  range.  However,  in  the  cases  where  the  stresses  corresponding  to 
minimum  load  become  negative,  the  stress  range  should  not  be  used.  Rather 
the  direct  input  of  maximum  and  minimum  stresses  should  be  used  and  care 
must  be  exercised  to  neglect  the  negative  stress  intensity  factors  corre¬ 
sponding  to  minimum  load  while  computing  the  ranges.  The  computed  values  of 
stress  intensity  factor  ranges  for  steel  and  aluminum  lugs  are  presented  in 
Figures  3-59  and  3-60,  respectively.  The  results  are  given  for  Rq/R^  ratios 
of  1.50,  2.25  and  3.00  and  far-field  stress  ratios  of  0.5  and  0.1  Since  this 
was  a  nonlinear  analysis,  the  far-field  stress  ratio  (ratio  of  minimum  to 
maximum  far-field  stress)  differs  from  the  crack-tip  stress  ratio  (ratio  of 
minimum  to  maximum  stress  intensity  factor).  Furthermore,  because  of  the 
nonlinearity,  no  effort  was  made  to  normalise  the  stress  intensity  factor 
ranges  in  Figures  3-59  and  3-60.  However,  the  crack  length  in  the  presenta¬ 
tion  has  been  normalized  by  the  net  ligament  of  each  lug. 

For  normalized  crack  lengths  less  than  about  0.35,  the  stress  Intensity 
factor  ranges  are  lower  for  an  Rq/ R^  ratio  of  1.50  and  higher  for  increasing 
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figure  Stress  Intensity  factor 


Ranges  for  Steel  Lugs 


R  /K  ratios.  For  normalized  crack  lengths  above  0.35,  the  trend  changes  to 
o  i 

higher  stress  intensity  factor  ranges  for  lower  kQ/Ri  ratios.  However, 

expressed  in  terms  of  unnormalized  physical  crack  length,  c,  the  stress 

intensity  factor  for  any  given  crack  length  is  always  higher  for  a  lower 

R  /R.  ratio, 
o  1 

To  assess  and  to  have  a  feel  for  the  effects  of  plastic  yielding  by 
comparing  with  linear  analysis,  the  stress  intensity  factor  ranges  are  nor¬ 
malized  by  ACTQ\/¥c,  where  a ao  is  fche  corresponding  far-field  stress  range. 
These  normalized  values  given  in  Table  3-13  have  limited  significance  in 
the  sense  that  these  values  are  applicable  to  the  specific  material  and  load¬ 
ing  condition.  The  values  can  be  directly  compared  with  the  normalized 
stress  intensity  factor  ranges  given  in  Table  3-3  (applicable  for  R=0.1  and 
R=0.5),  which  correspond  to  linear  analysis,  to  evaluate  the  effect  of  plas¬ 
ticity.  Table  3-3  can  be  assumed  to  be  applicable  for  both  materials, 
since  the  differences  were  less  than  3  percent. 

Comparison  of  these  two  tables  indicates  that  the  normalized  values  of  A  K 

for  the  elasto-plastic  analysis  are  lower  for  smaller  c/fR^-R^)  ratios  and  become 

higher  as  the  crack  grows  through  the  lug  ligament.  This  is  because  of  the 

stress  redistribution  in  the  lug  due  to  yielding.  In  the  yield  zone  and  its 

neighborhood,  the  stresses  will  be  lower  when  compared  with  the  elastic 

stresses.  Away  from  the  plastic  zone,  the  stresses  will  be  higher  than  elastic 

stresses,  to  satisfy  the  equilibrium.  The  above  behavior  is  true,  in  general; 

however,  for  an  R  /R.  ratio  of  1.50  and  a  stress  ratio  of  0.1,  the  nonlinear 
o  1 

values  are  consistently  lower  than  the  linear  values  for  almost  all 

c/(Rq-R^)  values.  The  reason  for  this  behavioral  change  can  be  explained 
as  follows.  For  this  particular  lug  configuration  and  stress  ratio,  there 
is  more  plastic  yielding  as  well  as  compressive  residual  stresses  near  the 
lug  hole  at  the  minimum  load.  These  negative  stresses  are  truncated  in  the 
stress  intensity  computation.  Also,  the  net  ligament  length  is  the  smallest. 

For  these  two  reasons,  i.e.,  truncation  of  negative  stresses  and  smaller 
Ligament  size,  the  normalized  values  of  stress  intensity  factor  range  never 
become  higher  than  those  of  the  linear  analysis.  Also,  note  that  due  to 
plastic  yielding  while  loading  and  subsequent  linear  unloading,  the  normalized 
values  of  k=  0.5  for  lower  c/fR^/R^)  ratios  are  greater  than  the  corresponding 
values  of  R =  ^ . 1  and  closer  to  linear  elastic  values  given  in  Table  3-3. 
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The  elasto-plastic  method  can  also  be  used  to  calculate  the  stress  inten¬ 
sity  factor  solutions  for  corner  cracks  emanating  from  attachment  lugs  that 
ai->  subjected  to  loadings  above  yield.  The  above  through-the-thicklesa  solu¬ 
tion  ran  bo  modified  with  front-free  surface,  curvature  correction,  etc.  for 
corner  crack  problems,  as  in  Equations  (34)  and  (35)  of  subsection  3. 
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SECTION  IV 


STRESS  INTENSITY  FACTORS  FOR  TAPERED  ATTACHMENT  LUCS 


Tapered  attachment  tugs  have  been  uced  frequently  in  aircraft  struc¬ 
tural  fittings  to  provide  strength  against  off-axis  loading.  To 
determine  the  critical  location  and  direction  where  a  fatigue  crack  may 
initiate  and  subsequently  grow,  it  is  necessary  to  determine  the  stress  dis¬ 
tribution  of  the  unflawed  lug  and  the  stress  intensity  factors  of  the  cracked 
lug.  This  section  describes  the  development  of  stress  distribution  and 
stress  intensity  factors  for  tapered  attachment  lugs  subjected  to  symmetric, 
off-axis  and  transverse  loadings. 

U  STRESS  ANALYSIS 

Because  of  the  complexity  of  the  lug  geometry  and  off-axis  loading, 
the  finite  element  method  was  used  to  determine  the  tangential  stress  dis¬ 
tributions  along  the  edge  of  the  hole  in  the  unflawed  lug.  This  analysis 
was  conducted  for  a  pin  load  applied  at  0°,  90°,  135°,  180°,  270°,  and  315° 
measured  in  the  clockwise  direction  from  the  axis  of  the  lug.  Figure  4-1 
depicts  the  geometry  and  typical  two-dimensional  finite  element  model  used 
in  the  stress  analysis  of  the  unflawed  lug.  The  angle  between  the  two  edge 
surfaces  of  the  tapered  head,  j3,  is  45  degrees.  The  finite  element  model 
shown  in  Figure  4-1  consists  of  429  nodes,  72  triangular  elements,  348  quad¬ 
rilateral  elements,  28  spring  elements,  and  a  total  of  850  degrees  of  free¬ 
dom.  A  concentrated  force  was  applied  at  the  center  of  the  pin  to  simulate 
pin  loading  and  was  reacted  at  the  base  of  the  lug.  The  analysis  was 
carried  out  to  determine  the  stresses  in  the  lug  and  the  pin-bearing  pres¬ 
sure  distributions  at  the  pin-lug  contact  area  for  the  six  loading  direc¬ 
tions  mentioned  above.  Three  outer-to-inner  radius  ratios  of  the  tapered 

head,  i.e.,  R  /K.  1.50,  2.25  and  3.00,  were  evaluated.  In  all  cases  of 

o  i 

the  analysis,  the  rigidity  oi  the  pin  was  assumed  to  be  three  times  the 
rigidity  of  the  lug,  which  simulates  an  aluminum  lug  loaded  by  a  steel  pin. 

For  the  convent iona l  finite-element  analysis  using  constant  strain 
elements,  the  values  o!  stress  and  strain  obtained  for  a  given  element 
were  assigned  to  the  centroid  location  of  that  element.  To  determine 
the  stress  at  the  edge  ot  the  hole,  the  stresses  at  the  centroids  of  a 
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Figure  4-1.  Finite  Element  Model  for  a  Tapered  Attachment  Lug 

Having  a  K  /R  Ratio  of  2.2^ 
o  i 


series  of  elements  located  in  the  same  radial  direction  were  used  to  extrap¬ 
olate  to  the  edge  location.  This  procedure  was  used  to  determine  the  tan¬ 
gential  stresses  along  the  edge  of  the  hole  of  tapered  lugs  having  Rq/R^ 
ratios  of  1.50,  2.25  and  3.00.  The  results  are  shown  in  Figures,  4-2,  4-3, 
and  4-4  for  pin  loadings  applied  in  three  principal  directions,  i.e.,  0°, 

-45°  (315°),  and  -90°  (270°),  respectively.  In  each  case,  the  tangential 
stresses  were  normalized  with  the  average  pin  bearing  pressure,  cr^,  which 
is  defined  as  P/(2R^B),  where  B  is  the  thickness  of  the  lug. 

As  can  be  seen  from  these  figures,  for  each  loading  direction,  there 
are  two  local  maximum  tangential  stresses  located  at  each  side  of  the  load¬ 
ing  direction.  The  locations  of  these  maximum  stresses  depend  upon  the 
loading  direction  and  Rq/ R ^  ratio.  For  a  pin  loading  applied  in  the  axial 
direction  of  the  lug,  as  shown  in  Figure  4-2,  the  maximum  stress  locations 
are  found  at  about  +85  to  +90°  away  from  the  loading  direction.  The  maximum 
stress  locations  do  not  change  significantly  with  the  change  of  the  Rq/R^ 
ratio.  When  the  pin  loading  is  applied  in  the  -45°  direction  (see  Figure 
4-3),  for  Rq/R^  =  1.50,  the  absolute  maximum  stress  occurs  at  about  65° 
measured  from  the  axis  of  the  lug  (or  110°  away  from  the  load  direction) 
with  the  other  local  maximum  stress  located  about  180°  from  the  first.  When 
the  Bj/Ki  ratio  increases,  the  locations  of  the  maximum  stresses  change 
only  slightly.  However,  the  absolute  maximum  stress  location  switches  from 
the  head  side  of  the  lug  to  the  base  side  of  the  lug.  When  the  pin  loading 
is  applied  in  the  direction  perpendicular  to  the  axis  of  the  lug  (see  Figure 

4-4) t  for  all  K  /R.  ratios,  the  maximum  stress  occurs  at  the  location  in 
o  i 

the  base  of  the  lug  at  a  0  value  of  about  200°  to  210°.  Figure  4-5  summa¬ 
rizes  the  locations  of  the  local  maximum  tangential  stresses  at  the  edge  of 
the  hole  for  each  loading  direction.  These  are  the  most  critical  locations, 
where  one  would  anticipate  that  a  fatigue  crack  would  initiate.  For  each 
lug  geometry  and  loading  condition.  Numbers  1  and  2  shown  in  Figure  4-5 
indicate  the  probable  order  of  crack  initiation.  They  arc  chosen  based  upon 
the  relative  magnitude  of  the  two  computed  local  maximum  stresses.  Stress 

distributions  along  the  x-axis  for  tapered  lugs  with  R  /R,  ranging  from  1.5  to 

o  t 

3.0  subjected  to  symmetric  loading  are  presented  in  Figure  4-6  and  Table  4-1. 

Larsson  [47,  48]  conducted  fatigue  testing  of  axially  and  transversely 

loaded  aluminum  lugs  having  a  R  /R.  ratio  of  2.2.  He  observed  and  tabu- 

o  t 

lated  the  locations  of  fretting  and  crack  initiation  for  three  different  lug 
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Figure  4-2.  Normalised  Tangential  Streaa  Distribution*  Along 
the  Edge  of  the  Hole  for  Tapered  Lugs  Subjected 
to  a  Pin  Loading  Applied  In  0®  Direction 
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Figure  4-4.  Normalised  Tangential  stress  Distr ibut ions  Along 
the  Edge  of  Che  Hole  for  Tapered  hug*  Subjected 
to  a  Fin  Loading  Applied  in  -<K>  Direction 


Figure  <*-h.  Unflawed  Stress  Distributions  at  Tapered  Attachment  hugs  Having  Various  K  /K  Ratios 


TABLE  4-1 NORMALIZED  UNFLAWED  STRESS  DISTRIBUTION  ALONG  x-AXIS  POR  TAPERED  ATTACHMENT  LUGS 


♦Extrapolated  Values 


configurations.  His  results  of  crack  initiation  locations  are  summarized  in 
Figures  4-7  and  4-8  for  a  pin  loading  applied  in  the  0°  and  90°  directions, 
respectively.  As  seen  from  these  figures,  the  current  predicted  fatigue 
critical  locations  as  well  as  the  possible  sequence  of  crack  Initia¬ 
tion  agree  well  with  these  experimental  data.  These  predicted  critical 
locations  will  be  used  in  the  modeling  of  the  cracked  lugs  in  the  subsequent 

iracture  analysis. 


The  fatigue  critical  area  of  an  attachment  lug  for  a  given  loading  di¬ 
rection  is  not  necessarily  subjected  to  compression  when  the  loading  direc¬ 
tion  is  reversed.  Thus,  for  some  load  orientations,  the  reversed  fatigue 
loading  might  have  a  significant  effect  on  crack  growth  behavior.  Figures 
4-9  through  4-11  show  the  tangential  stress  distributions  along  the  edge  of 
the  hole  for  a  pin  loading  applied  in  the  reversed  direction  of  the  three 
primary  load  orientations  presented  in  Figures  4-2  through  4-4, 
respectively.  As  seen  from  :  he  figures,  in  most  cases,  the  reversed 
loading  stretches  the  critLcal  area  in  tension,  but  the  magnitude  is 
reduced.  Note  that  the  result  In  F;  >ure  4-11  is  essentially  the  same 
as  that  in  Figure  4-4  except  for  the  definition  of  the  angle  0  .  Plots  of 
the  stress  concentration  factors  at  the  edge  of  the  hole  in  logarithmic 
scales,  as  shown  in  Figure  4-12  for  symmetrically  loaded  tapered  attachment 
lugs, reveal  that  a  simple  empirical  formula  can  be  derived  and  the  relation¬ 


ship  is  given  by  the  equation 
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where  is  the  taper  angle  of  the  attachment  lug  in  degrees.  Note  that 

tor  straight  lugs  (/I  0°)  the  above  equation  reduces  to  that  given  in 

Figure  3~*>. 


The  stress  concentration  factors  were  obtained  by  normalizing  the  peak 
tangential  stresses  wich  the  average  pin  bearing  pressure,  (7  .  The  equa¬ 

tion  of  the  logarithmic  straight  lines  shown  in  Figure  4-1?  compares  with 
the  finite  element  solution  within  O.H  percent  for  tapered  lugs.  The  values 
computed  by  the  above  equation  and  the  finite  element  method  are  also  listed 
in  tabular  form  in  Figure  4-1?  for  comparison.  This  simple  empirical  equa¬ 
tion  may  be  used  for  interpolating  for  taper  angles  less  than  4'>°  ot  extrap¬ 
olating  tor  K  /H,  values  outside  the  range  o!  l.y  and  3.0.  The  values  of 
o  i 
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the  stress  concentration  factors  presented  in  Figure  4-12  correspond  to  the 
0  ~  _t-90°  locations.  The  actual  maximum  tangential  stress  location  is  be¬ 
tween  +85°  and  j-90°.  From  a  knowledge  of  the  stress  distributions,  it  is 
anticipated  that  there  will  be  no  appreciable  difference  in  the  solutions. 
Also,  in  the  fracture  analysis  of  tapered  lugs  subjected  to  a  pin  loading 
in  the  axial  direction,  the  crack  surface  will  be  assumed  to  be  in  the  plane 
perpendicular  to  the  loading  direction.  However,  for  off-axis  loadings  of 
+45°  and  ^-90°,  the  crack  surface  will  be  assumed  to  be  in  the  critical  loca¬ 
tion  as  predicted  by  the  stress  analysis. 

The  computed  pin  bearing  pressure  distributions  along  the  contact  sur¬ 
faces  with  no  cracks  present  are  presented  in  Figures  4-13  through  4-15  for 
pin  loadings  applied  in  the  three  primary  directions,  i.e.,  0°,  -45°,  and 
90°,  respectively.  Similar  results  obtained  lor  the  reversed  loadings  of 
180°  and  135°  are  presented  in  Figures  4-16  and  4-17,  respectively.  For  the 
reversed  loading  of  90°,  the  contact  pressure  can  be  obtained  from  Figure 
4-15  by  redefining  the  angle  ll  as  previously  discussed. 

2._  STRESS  IN’TKMSTTY  FACTOR  ANALYSIS 

based  on  the  stress  analysis  of  unflawed  tapered  attachment  lugs,  the 
most  critical  locations  were  selected  for  the  fracture  analysis  to  obtain 
the  stress  intensity  factors  for  various  crack  lengths.  The  special  high- 
order  crack  tip  singularity  element  [8]  was  again  used  in  the  present  analy¬ 
sis  to  calculate  the  stress  intensity  factors.  In  the  analysis,  it  is 
assumed  that,  tor  a  given  tapered  lug  subjected  to  a  specific  direction  of 
pin  loading,  the  crack  will  initiate  from  the  maximum  tangential  stress 
location  at  the  hole  and  propagate  radially.  Figure  4-18  shows  a  typical 
finite  element  model  which  was  used  for  a  single  through-the-th ickness  crack 
emanating  from  a  tapered  attachment  lug  subjected  to  a  pin-loading  applied 
in  the  0°  and  180°  directions.  The  crack  surface  is  assumed  to  be  in  the 
plane  perpendicular  to  the  loading  direction.  The  finite  element  breakdown 
consists  of  366  nodes,  150  triangular  elements,  228  quadrilateral  elements, 
32  spring  elements,  1  crack-tip  element  ami  a  total  of  724  degrees  of  free¬ 
dom.  The  computed  normalised  opening  mode  stress  intensity  factors,  as  a 
function  of  normalized  crack  length,  for  single  through-the-thickness  cracks 
In  tapered  attachment  lugs  having  !<  /k.  ratios  ranging  1  rom  1.5 

O  1 

Co  1.0  are  presented  in  Figure  4-14  and  Table  4-2*  The  normalized  stress 
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Figure  4-13.  Pin-Bearing  Pressure  Distributions  Along  the 

Contact  Surface  of  Tapered  Lugs  Subjected  to  a 
f  Pin  Loading  Applied  in  0°  Direction 
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Figure  4-15.  Pin-Bearing  Pressure  Distributions  Along  the 

Contact  Surface  of  Tapered  Lugs  Subjected  to  a 
Pin  Loading  Applied  in  -90°  Direction 
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.  Pin-Hearing  Pressure  Distributions  Along  the 
Contact  Surface  of  Tapered  hugs  Subjected  to  a 
Pin  Loading  Applied  in  180°  Direction 
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Figure  4-17.  Pin-3earing  Pressure  Distributions  Along  the 

Contact  Surface  of  Tapered  Lugs  Subjected  to  a 
Pin  Loading  Applied  in  135°  Direction 
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Figure  4-19,  normalized  Stress  Intensity  Factors  for  Single  Through— the-Thickness  Cracks 
Emanating  from  Tapered  Attachment  hugs  Subjected  to  a  Pin  Loading  Applied 
in  0  Loading  Direction 


TABLE  4-2.  NORMALIZED  STRESS  INTENSITY  FACTORS  FOR  SINGLE  THROUGH— THE— THICKNESS  CRACKS 
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intensity  factor  values  at  the  lug  hole,  i.e.,  c/R^  =  0,  were  obtained  by 
multiplying  the  concentration  factors  determined  from  the  unflawed  stress 
analysis  by  1.12,  which  was  derived  by  Gross  et  al  [ 2 1 ]  for  a  straight  edge 
crack  in  a  large  plate  subjected  to  remote  tension.  The  trends  of  the  above 
tapered  lug  solutions  are  very  similar  to  those  obtained  for  straight  lugs. 
Similar  solutions  of  normalized  stress  intensity  factors  as  a  function  of 
normalized  crack  length  for  tapered  attachment  lugs  subjected  to  pin  load¬ 
ings  applied  in  the  180°  direction  are  presented  in  Figure  4-20. 

The  analysis  was  then  extended  to  tapered  attachment  lugs  subjected  to 
off-axis  loadings  of  -45°,  135°,  -90°  and  90°.  A  R  /R^  ratio  of  2.25  was 
considered  for  the  off-axis  loading  fracture  analysis.  Figure  4-21  shows 
the  finite  element  model  used  for  analyzing  a  single  through-the-thickness 
crack  emanating  from  a  tapered  lug  subjected  to  a  pin  loading  applied  in 
-45°  and  its  reversed  (135°)  direction.  The  model  contains  401  nodes,  224 
triangular  elements,  218  quadrilateral  elements,  34  spring  elements,  1  crack- 
tip  element  and  a  total  of  794  degrees  of  freedom.  Note  that  although 
there  are  two  critical  locations  modeled,  58°  and  227°  measured  from 
the  axis  of  the  lug,  only  one  crack  was  analyzed  at  a  time.  The  computed 
normalized  stress  intensity  factors  are  shown  in  Figure  4-22  as  a  function  of 
the  normalized  crack  length.  This  figure  shows  that  when  the  crack  length 
is  small,  say  c/ Ri  <  0.15,  the  stress  intensity  factor  for  a  crack  located 
closer  to  the  base  of  the  lug  is  higher  than  the  one  located  at  the  head  side 
of  the  lug.  When  the  crack  length  increases  (c/R^  >  0.15),  the  stress  inten¬ 
sity  factor  for  a  crack  located  at  the  head  side  becomes  larger  than  the  one 
located  at  the  opposite  side  of  the  hole.  The  difference  between  the  two 
computed  K-values  increases  as  the  crack  length  increases.  For  the  case  when 
the  direction  of  the  applied  pin  loading  is  reversed,  the  stress  intensity 
factors  are  also  computed  and  included  in  the  figure  for  a  crack  located  at 
the  head  side  of  the  lug.  The  surfaces  of  a  crack  located  closer  to  the  base 
of  the  lug  are  completely  closed  during  the  reversed  loading.  The  computed 
K-values  corresponding  to  the  reversed  loading  are  much  smaller  than  the 
corresponding  ones  obtained  under  the  primary  tensile  loading.  In  the  above 
analysis,  only  one  crack  was  assumed  to  exist  at  a  time.  However,  in  reality 
both  cracks  may  exist  and  grow  at  the  same  time.  In  such  a  case,  the  influ¬ 
ence  of  one  crack  on  the  stress  intensity  factor  of  the  other  crack  may  be 
significant.  This  cart  be  accounted  for  by  developing  a  matrix  of  stress 
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Figure  4-2 1. 


Finite  K  lenient  Model  for  a  Cracked  Tapered  Lug 
Subjected  to  a  Pin  Loading  Applied  in 
-45  and  135°  Directions 
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NORMALIZED  STRESS  INTENSITY  FACTOR,  K/(o 
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Figure  4-22.  Normalized  Stress  Intensity  Factors  tor  Single 
Through-the-Thlckness  Crack’  Emanating  from  a 
Tapered  Lug  Subjected  to  a  Pin  Loading  Applied 
in  -45  and  its  Reversed  Directions 


intensity  factor  solutions  for  various  crack  lengths  of  both  cracks.  This 
matrix  of  solutions  may  then  be  used  to  accurately  estimate  the  stress  in¬ 
tensity  factor  for  either  crack.  However,  such  an  effort  is  not  made  in  the 
present  analysis. 

A  similar  finite  element  model  to  the  one  shown  in  Figure  4-21  was  used 
for  a  tapered  lug  subjected  to  a  pin-loading  applied  in  the  +-90°  directions. 

Two  critical  locations,  43°  and  205°  measured  from  the  axis  of  the  lug,  were 
determined  from  the  results  of  the  stress  analysis  of  an  unflawed  tapered 
lug  subjected  to  a  pin  loading  applied  in  the  -90°  direction.  Similar  to 

the  -45°  loading  case,  only  one  crack  was  analyzed  at  a  time,  though  two 

critical  locations  were  modeled.  The  computed  normalized  stress  intensity 
factors  are  shown  in  Figure  4-23  as  a  function  of  the  normalized  crack  length. 
The  figure  shows  that  when  the  crack  length  is  smaller  than  0.85  , 

the  stress  intensity  factor  for  a  crack  located  closer  to  the  base  of  the  lug 
is  higher  than  the  one  located  at  the  head  side  of  the  lug.  When  the  crack 
length  increases  (c/K^  >  0.85),  the  stress  intensity  factor  for  a  crack 
located  at  the  head  side  becomes  larger  than  the  one  located  at  the  opposite 
side  of  the  hole.  For  the  case  when  the  direction  of  the  applied  pin  loading 
is  reversed,  the  stress  intensity  factors  are  also  computed  and  included  in 

the  l  igure  tor  a  crack  located  at  the  head  side  of  the  lug.  The  surfaces  of 

a  crack  located  closer  to  the  base  of  the  lug  are  completely  closed  during 
the  reversed  loading.  The  K-values  corresponding  to  the  reversed  loading  are 
smaller  than  the  corresponding  ones  obtained  under  primary  tensile  loading. 
However,  these  magnitudes  are  significantly  larger  than  the  corresponding 
ones  obtained  tor  a  case  where  the  pin  loading  was  applied  in  the  135°  direc¬ 
tion. 

Based  on  the  current  unflawed  stress  and  stress  Intensity  factor  analysis, 

one  stay  conclude  that,  among  the  principal  pin  loading  directions  of  0°,  -45° 

and -90*' ,  the  loading  of  -45°  is  the  most  severe  loading  case  from  the  point  of 

view  of  fatigue  crack  initiation  and  fatigue  crack  propagation.  Also,  when  a 

tapered  I  ue  havinv  a  R  /R.  ratio  of  2.25  is  subjected  to  cyclic  fatigue  pin 

o  1 

loading  in  -45'  and  -91)  directions,  a  crack  will  probably  initiate  at  the 
critical  location  closer  t o  the  base  ot  the  lug  first.  When  this  crack 
propagates,  a  second  crack  will  initiate  at  the  head  side  of  the  lug.  Eventu¬ 
ally,  the  g  t\>wt  h  of  the  base-side  crack  will  slow  down  while  the  growth,  rate  of 
the  head-side  crack  will  increase  and  exceed  that  of  the  base-side  crack.  Tile 
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NORMALIZED  STRESS  INTENSITY  FACTOR,  K/(crbr  JTc) 


NORMALIZED  CRACK  LENGTH,  c/R. 


Figure  A-?).  Normalized  Stress  Intensity  Factors  for  Single 
Through-the-Thickness  Cracks  Emanating  from  a 
Tapered  t.ug  Subjected  to  a  Pin  Loading  Applied 
in  -90  and  its  Reversed  Directions 
<K  /K,  .  2.25) 
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lug  will  finally  fail  from  the  head  side.  Failure  inodes  experimentally 
studied  by  Larsson  [46]  for  tapered  lugs  loaded  in  the  -90°  direction  are 
very  similar  to  those  discussed  above. 


SECTION  V 


CRACK  GROWTH  ANALYSIS  METHOD 


Linear  elastic  fracture  mechanics  methodology  is  discussed  in  this 
section  to  predict  crack  growth  in  damaged  attachment  lugs.  An  accurate  and 
efficient  prediction  of  crack  growth  from  an  initial  crack  length  to  a  criti¬ 
cal  or  final  crack  length, when  subjected  to  a  prescribed  load  history,  de¬ 
pends  on  a  damage  accumulation  package  containing  the  following  four  basic 
elements: 

o  Applied  load  spectrum 

o  Stress  intensity  factor  solution 

o  Baseline  crack  growth  rate  data  or  crack  growth  rate  equation 

o  Spectrum  load  interaction  model. 

The  applied  load  spectrum  is  the  loading  sequence  or  history  the  attach¬ 
ment  lug  is  subjected  to  for  which  the  crack  growth  prediction  is  to  be  made, 
and  is  therefore  basically  an  input  for  the  problem. 

ft  is  generally  accepted  that  the  stress  intensity  factor  controls  the 
rate  of  propagation  of  a  fatigue  crack  in  structural  components.  Any  cracked 
structure  can  be  expected  to  respond  in  a  predictable  manner  to  a  given 
applied  load  spectrum  if  its  associated  stress  intensity  factor  solution  is 
available.  The  development  of  such  stress  intensity  factor  solutions  was  pre¬ 
sented  in  the  previous  two  sections  for  straight  and  tapered  attachment  lugs 
for  various  configuration,  design  and  loading  complexities. 

The  response  of  the  crack  is  usually  given  in  the  form  of  constant 
amplitude  fatigue  crack  propagation  rate,  da/dN,  as  a  function  of  stress 
intensity  factor  range,  ZiK.  These  baseline  material  proDertv  data  are 
obtained  exnpr imental Iv  and  can  Ke  exnressed  in  either  tabular  or  equation 
form,  i .e . , 


~  i  (AK,  . )  (4b) 

uN 

Paris’  [4d],  Forman's  t  *>0  ]  and  Walker's  [m]  eq:ations  are  some  of  the 
typical  representative  crack  growth-rate  equations  which  will  be  discussed 
in  this  section. 


The  normal  crack  growth  rate  under  constant  amplitude  loading  changes 
if  the  load  cycle  is  preceded  by  a  different  amplitude  load  cycle.  As  ex¬ 
amples,  a  tensile  overload  causes  permanent  plastic  deformation  at  the  crack 
tip  which  in  turn  delays  the  crack  growth  at  subsequent  low-load  cycles, 
while  the  compressive  underload  may  accelerate  the  normal  crack  growth. 

To  account  for  such  spectrum  load-interaction  effects,  several  crack- 
growth  retardation  models  have  been  proposed.  Models  such  as  Wheeler  [52], 
Wilienborg  [53],  Generalized  Wiilenborg  [54]  and  Hsu  [55]  will  be  discussed 
in  this  section. 


CRACK  GROWTH  RATE  EQUATIONS 


A  large  and  growing  volume  of  crack  growth  data  are  available.  Most  of  the 
data  are  expressed  in  terms  of  AK  versus  da/dN.  In  this  relationship, 
there  is  a  threshold  value  of  K  ,  say  K^,  heiow  which  a  flawed  structure 
can  be  cycled  without  measurable  crack  extension.  At  the  other  extreme  is 
Kc  ,  called  fracture  toughness,  a  value  of  K  at  which  a  flaw  will  propa¬ 
gate  unstably.  Therefore,  fatigue  crack  growth  must  occur  in  the  range 

K  ,  to  K  .  All  crack  growth  rate  data  available  for  establishing  the 
the 

da/dN  versus  AK  relationship  fall  within  this  range.  Normally,  data 
show  some  experimental  scatter,  and  the  least-square  mean  line  representa¬ 
tion  of  all  test  data  is  used  to  establish  the  AK  versus  da/dN  relation¬ 
ship. 


Several  methods  of  using  the  mean  data  to  describe  crack  growth  rates 
have  been  proposed.  Of  these  methods,  three  well  established 

crack  propagation  equations  are  used  most.  Thev  are  proposed 
by  Paris  et  al  [49],  Forman  et  al  [50  ],  and  Walker  [51  ].  Paris*  equation  is 
limited  to  cases  where  a  constant  stress  ratio,  K,  is  applied.  Forman  and 
Walker's  equations  are  more  general  and  applicable  to  cases  where  variable 
loads  exist  in  the  load  spectrum. 


Accuracy  of  the  predicted  crack  growth  life,  using  any  crack  growth 
equation,  depends  upon  accurate  values  of  da/dN  versus  AK  and  correct 
stress  intensity  factor  values  for  a  particular  geometry. 


1 . 1  Paris'  Equation 

The  crack  growth  rate  equation  proposed  by  Paris,  et  al,  is  of  the 
tallowing  form: 


(47) 


where: 


n 

da/dN  -  c  (AK)  ^ 
P 


da/dN  =  rate  of  crack  growth  (inch/cycle) 

c  and  n  are  constants  for  a  particular  stress  ratio 
P  P 


AK  -  stress  intensity  factor  range  (K  -  K  ,  ). 

max  min 

Paris'  equation  has  been  an  effective  regression  equation  for  crack 

propagation  rate  in  unreinforced  flat  specimens  under  constant  amplitude 

loading.  However,  it  is  not  successful  in  pooling  such  data  for  more  than 

one  value  of  stress  ratio  K,  (n  /n  ). 

mm  max 


1.2  Forman's  Equation 

In  an  effort  to  extend  Paris'  equation  to  cover  various  stress  ratios, 
and  to  take  into  account  the  instability  of  the  crack  growth  when  the  stress 
intensity  factor  approaches  its  critical  value  for  a  given  material,  Forman, 
et  al, modified  Paris'  equation  as  follows: 


where: 


da/dN 


cf(AK) 

(  i-uk  -  AK 
c 


(4H) 


c.  and  n,  -  material  constants 
f  f 

K  critical  stress  intensity  factor 
c  1 

because  of  crack-tip  blunting  which  is  produced  by  cyclic  strains  in 

the  material  ahead  ot  the  advancing  crack,  it  is  important  to  select  the 

appropriate  K  value.  The  choice  of  a  K  value  larger  than  the  fracture 

^  c 

toughness  ot  the  material  normally  will  give  a  better  result  in  establishing 
K  versus  da/dN  relationships.  It  the  available  data  includes  more  than 
one  k-ratio,  should  he  used  as  a  curve  fitting  parameter  to  force  the 

best  fit  £‘>r  all  available  data. 


1.3  Walkef's  Equat ion 

Roberts  and  Erdogan  [36]  proposed  that  stress  ratio  effects  in  the  power 
law  region  of  the  crack  growth  rate  curve  could  be  described  using 

P  P 

da/dN  »  c  (K  1  .  AK  2)  (49) 

max 

Walker  [51]  reformulated  the  equation  so  that  the  stress  ratio  would  be 
explicit  and  it  had  the  form 


da/dN  -  c 


where  c,  m  and  p  are  material  constants. 


(50) 


_l_.4_  Determination  of  Constants  in  Paris'  and  Forman’s  Equations 

Taking  logarithms  of  both  sides  of  Paris'  and  Forman's  equations  and  re¬ 
arranging  terras,  they  become,  respectively, 

log  (da/dN)  a  log  c^  ♦  n^  log  (AK)  (31) 

and 

log  j  [a  -  R)  Kc  -  Akj  da/dN  ^  log  cf  ♦  nf  log  (AK)  (52) 
For  any  two  given  coordinate  points,  say  AK  ,  da/dNi  and  Ak  , 

^  1  4  i  +  1 

la/ dN i ^ j  which  represent  a  segment  of  the  growth  rate  curve,  one  can  solve 
the  two  simultaneous  equations  for  the  unkovns  c  and  n. 

2.  SPECTRUM  LOAD  INTfcw..  .TIQN  MODELS 

The  majority  of  crack  growth  studies  have  been  carried  out  for  constant 
stress  amplitude  because  of  the  simplicity  in  testing  and  because  the  data 
can  be  presented  in  a  straightforward  manner.  In  actual  practice,  engineer¬ 
ing  components  are  subjected  to  loads  which  are  often  irregular  and  varv  in 
a  random  manner.  Fluctuating  loads  due  to  atmospheric  turbulence  on  aircraft 
are  typical  examples  of  random  loadings.  The  difference  between  the  conven¬ 
tional  fatigue  experiments  and  the  actual  loading  of  the  components  has  led  t 
the  development  of  various  crack  analysis  approaches.  As  discussed  ir  a  pre¬ 
vious  section,  normal  crack  growth  rate  under  constant  amplitude  loading 
changes  if  the  load  application  is  preceded  by  a  loading  cycle  of  a  different 
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amplitude.  The  tensile  overload  causes  permanent  plastic  deformation  at  the 
crack  tip  which  in  turn  delays  the  crack  grcwth  at  subsequent  low  load  cycles, 
while  the  compressive  overload  may  accelerate  the  normal  crack  growth.  The 
importance  of  delay  in  the  rate  of  fatigue  crack  growth,  as  produced  by  ten¬ 
sile  overloads,  on  the  accurate  prediction  of  fatigue  lives  of  structures  has 
been  recognized,  and  quite  a  few  investigations  have  been  stimulated  in  this 
area. 


Several  models  have  been  suggested  to  account  for  the  effects  of  delay 
on  the  prediction  ol  fatigue  crack  growth,.  Those  models  use  the  plastic- 
zone  size  (for  either  plane  stress  or  plane  strain)  associated  with  the 
applied  load  levels  to  characterize  the  load-interaction  effect.  These 
models  assume  that  if  the  size  of  the  plastic  zone,  r  ,  developed  due  to 
application  of  the  current  load  cycle  at  crack  length  a  extends  to  or  past 
the  extremities  of  a  previously  developed  load-interaction  zone,  a^,  (i.e., 
(a  +  r  )  >  a  ),  there  will  be  no  load  interaction  and  the  growth  rate  asso- 
ciated  with  the  current  load  cycle  is  the  same  as  the  one  generated  under 
constant-amplitude  loading.  Conversely,  the  crack-growth  rate  would  be  re¬ 
duced  (retarded)  if  r  <  (a  -  a).  Some  of  the  most  often  referenced 

P  P 

retardation  models  are  described  below. 


2.1  Wheeler  Model 

The  first  retardation  model  proposed  was  that  of  Wheeler  [52]  who  sug¬ 
gested  the  use  of  a  reduction  factor  on  the  constant  amplitude  crack  growth 
rate: 


da 

dN  Cp 


f  <AK) 


where  c  is  a  retardation  parameter  defined  as 


for  a  -  a  >  K 

P  y 


c  -  t 
P 


tor  a  -  a  <  K 

p  y 


(53) 
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and  m  =  shaping  exponent 

Ry  -  extent  of  the  current  yield  zone 

a^  -  a  =  the  distance  from  crack  tip  to  elastic-plastic  interface 

This  model  requires  previous  spectrum  growth  data  to  derive  an  empirical 
"shaping  exponent,"  m.  Although  it  is  a  substantial  improvement  over  the 
linear  cumulative  damage  rule,  this  model  is  more  of  a  data  fitting  technique 
than  it  is  a  predictive  technique.  Besides,  the  exponent  in  Wheeler’s  re¬ 
tardation  parameter  is  dependent  not  only  upon  the  material,  but  also  upon 
the  manner  in  which  spectrum  loads  were  applied.  For  any  design  spectrum 
different  from  the  ones  used  to  generate  the  test  data,  the  exponent  m 
may  be  very  different. 


2.2  Willenborp.  Model 

As  mentioned  earlier,  the  rate  of  propagation  of  a  fatigue  crack  is 
controlled  by  the  stress-intensity  factor  at  the  crack  tip.  Therefore,  the 
magnitude  of  the  stress-intensity  factor  is  a  good  indicator  of  the  extent 
of  crack  tip  deformation.  Since  the  cyclic  crack  tip  deformation  is  reduced 
doe  to  prior  overloads,  the  model  of  Willenborg  et  al  [53],  which  uses  an 
"effective  stress"  concept  to  reduce  the  applied  stress,  and  hence  the  crack 
tip  stress-intensity  factor,  seems  more  appealing.  The  model  is  described 
as 

S  '  1  [<AK)«ff  <K~,)eff]  <54) 

where 
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red  ap  max 
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(55) 


Jap 


a 

c 


and 


a  =  the  crack  length  at  the  beginning  of  the  load  cycle 

after  the  overload 

a  -  a  =  current  plastic  zone  caused  by  the  overload 
p  c  K  3 

F  =  tensile  yield  strength  of  the  material 

0  =  total  geometric  correction  factor 


This  model  can  be  used  to  predict  the  fatigue  crack  growth  under  spec¬ 
trum  Loading  without  the  assistance  of  empirical  factors  or  data.  Although  the 
model  gives  fairly  ,;ood  predictions  for  crack  growth  under  moderate  spectra, 
it  tends  not  to  give  good  correlations  for  spectrum  loading  with  high  over¬ 
loads,  especially  overloads  exceeding  150  percent  of  the  maximum  of  the  low 
load. 


2.1  Generalized  Willenborg  Model 

According  to  the  Willenborg  et  al  model,  if  the  overload  ratio, 

K  greater  than  or  equal  to  2.0,  a  zero  effective  r tress-intensitv 

UtcIX  lu-i  X 

factor  wLll  be  predicted  and  the  crack  will  stop  growing.  However,  experi¬ 
mental  data  contradicts  such  a  prediction.  Gallagher  [54]  modified  the 
Willenborg  equation  for  eifectLve  stress  Intensity  factors  as  follows: 
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or  zero,  whichever  is  larger. 
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The  factor  0  is  defined  as 
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t  -  (K 


max 


) 

Th 

-  1 


/K 

max 


(57) 


in  which  (K  )  is  the  maximum  threshold  stress  intensity  factor,  and 

max  Th  OL 

is  the  "overload  shut-off  ratio."  When  K  /K  >  crack  arrest 

SO  max  max  SO 

occurs.  Note  that  the  original  Willenborg  model  is  recovered  by  setting 

<t>  =  1. 

2.4  Hsu  Model 

The  Hsu  model  [55]  developed  at  tne  Lockheed-Georgia  Company  utilizes  an 
effective  stress  and  closure  concept.  It  assumes  that  the  stress  singularity 
does  not  exist  if  the  crack  surface  is  closed  and  that  the  crack  propagates  only 
during  that  portion  of  the  load  cycle  in  which  the  crack  surface  is  fully 
open.  Let  <jq  be  the  crack-opening  stress,  i.e.,  the  corresponding  far- 
fietd  stress  at  the  onset  of  crack  opening.  Then  the  effective  stress  ranee 
of  the  load  cycle  during  crack  propagation  can  be  defined  as 


ACT 


eff 


°max  °o 


(58) 


When  the  opening  stress  is  less  than  the  minimum  stress  of  the  applied 
load  cycle,  the  effective  stress  range  is  given  by 


ACT  re  =  q  -  <7  .  I  54  ) 

w  eff  max  rain 

If  the  maximum  stress  of  the  applied  load  cycle  is  less  than  the  crack 
opening  stress  due  to  prior  loads,  the  crack  surface  will  be  fully  closed. 
Hence,  the  fatigue  crack  will  not  propagate. 


However,  experimental  evidence  indicated  that  below  the  closure  K, 

strain  concentration  at  the  vicinity  of  the  crack  tip  still  exists.  Since 

fatigue  damage  is  normally  related  to  the  cyclic  strain  range,  the  effective 

K  ,  is  likely  ro  be  somewhat  lower  than  the  level  at  the  onsest  of  closure 
min 

and  higher  than  K  ^  under  steady  state  (constant  amplitude)  conditions. 
The  effective  stress  range  and  the  effective  load  ratio  can  then  be  re¬ 
written  as 
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A(Teff  ~  °"max  “  ^min*  ,c  (6°) 

eff 

Raff  =  (°min)  ff/(Tmax  <61> 

eff 

respectively,  where  cr  .  <  (G  .  )  s  rr  • 

nun  -  min  rr  S  v0 
eff 

For  any  given  load  cycle,  the  effective  stress-intensity  factor  can 
then  be  calculated  from  the  equation 

AK£ff  =Aaeff  .VfiT  .  ]3t  (62) 

The  crack  growth  rate  associated  with  this  applied  load  cycle  is  com¬ 
puted  from  the  growth  rate  equation 
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(63) 


For  a  given  loading  spectrum,  if  one  can  determine  the  effective  mini¬ 
mum  stress  corresponding  to  each  load  cycle,  the  fatigue  crack  propagation 
life  can  be  predicted. 


3.  CKACK  GROWTH  ANALYSIS  PROGRAM 

A  computer  program  has  been  developed  using  the  state-of-the-art 
methodologies  including  the  stress  intensity  factors  developed  under 
this  program  for  the  prediction  of  fatigue  crack  growth  behaviors  of  single 
through-the-thickness  cracks  and  single  corner  cracks  at  attachment  lugs 
under  cyclic  loading.  All  the  data  presented  in  Sections  III  and  IV,  such  as 
stress  analysis,  fracture  analysis.  Green’s  function,  interference-fit 
hushing  analysis,  etc.,  have  been  embedded  in  the  computer  program  to  make 
it  as  automatic  and  simple  as  possible.  This  computer  program  contains  the 
three  crack  growth  rate  equations  discussed  above.  There  are  five  different 
options  to  input  each  Individual  mission  load  profile  from  which  a  mission 
mix  spectrum  can  be  generated.  The  program  predicts  the  crack  growth  using  a 
block-by-block  integration  technique.  Crack  growth  may  be  analyzed  with  or 
without  load- interact  ion  using  any  of  rhe  models  described  above. 


For  through-the-thickness  cracks,  either  the  compounding  solution  or  the 
Green's  function  solution  can  be  used  in  the  prediction.  In  predicting  the 
growth  behavior  of  a  single  corner  crack,  the  crack  may  be  analyzed  by 
either  the  one-parameter  (i.e.,  constant  a/c  ratio)  or  two-parameter  method. 

For  one-parameter  analysis,  the  prediction  is  straightforward  and  is 
similar  to  through-the-thickness  crack  prediction.  For  two-parameter 
analysis,  it  is  assumed  that  for  a  given  number  of  applied  load  cycles,  the 
extension  of  the  quarter  elliptical  crack  border  is  controlled  by  the  stress 
intensity  factors  at  the  intersection  of  the  crack  periphery  at  the  hole 
wall  and  the  lug  surface,  i.e.,  K  and  K  ,  respectively.  In  general,  the 

A  v* 

stress  intensity  factors  at  these  two  locations  are  different,  resulting  in 
different  crack  growth  rates.  Therefore,  the  new  flaw  shape  aspect  ratio  after 
each  crack  growth  increment  will  be  different  from  the  preceding  one.  The  new 
crack  aspect  ratio  is  computed  using  the  new  crack  lengths  on  both  the  hole  wall 
and  lug  surface.  The  process  will  be  repeated  until  the  crack  length  along 
the  hole  wall  is  equal  to  the  lug  thickness.  At  that  time  the  transitional 
crack  growth  criteria  as  discussed  in  Section  III  are  used  until  the  crack 
has  achieved  a  uniform  format.  After  that,  if  the  failure  has  not  occurred, 
a  one-dimensional  through-the-thickness  crack  analysis  is  used  to  contin¬ 
uously  predict  the  subsequent  crack  growth  life.  The  analysis  is  considered 
to  be  complete  when  fracture  occurs  or  when  the  desired  final  crack  length 
or  the  maximum  usage  time  is  reached. 

A  detailed  description  of  the  computer  program  along  with  the  input 
instructions  is  provided  in  Ref erence  [57} # 


SECTION  VI 


SUMMARY  OF  RESULTS 


Analytical  methods  have  been  presented  to  predict  both  fatigue  crack 
growth  and  residual  strength  of  cracked  attachment  lugs.  Each  crack  growth 
analysis  includes  the  following  elements: 

o  Stress  intensity  factor  solution 
o  Baseline  crack  growth  rate  relationship 
o  Applied  load  sequence 
o  Spectrum  load  interaction  model. 

Of  these,  the  emphasis  in  this  report  has  been  upon  the  calculation 
of  stress  intensity  factors,  covered  in  Sections  III  and  IV.  These  sections 
have  been  followed  by  a  discussion  in  Section  V  of  the  alternative  constant 
amplitude  fatigue  crack  growth  rate  relationships  and  spectrum  load  inter- 
ac  t ion  mode  1 s . 

The  following  summary  paragraphs  are  intended  to  provide  an  overview 
and  to  tie  together  the  content  of  this  report  and  aid  the  reader  in  using 
the  analytical  methods  and  results  that  have  been  presented. 

1 •  NUMMARY  OF  STRESS  INTENSITY  FACTORS  FOR  STRAIGHT  LUCS  (SECTION  III) 

Section  III  has  covered  the  calculation  of  stress  intensity  factors 
tor  straight  attacluaent  lugs  subjected  to  axial  loading.  Several  alternative 
methods  have  been  discussed,  including  the  simple  compounding,  two-dimensional 
cracked  finite  element,  weighting  function,  and  three-dimensional  cracked 
finite  element  method.  Parameters  and  complexities  covered  in  the  stress 
intensity  factor  solutions  presented  in  this  section  are  outer-to-iimer  radius 
ratio  ( l . '0  to  1.0),  eraek  geometry  (single  corner  eracn,  through-the-thick- 
ness  crack,  and  the  intermediate  t  rails  i  t  i  on ),  e  r  ack  length  (measured  on  lug 
face  and  along  bore  of  hole),  change  in  distribution  of  pin  bearing  pressure 
due  to  eraek  length  change,  ratio  of  pin  modulus  to  lug  modulus  (1.0  or  1.0), 
interference-fit  bushings,  and  e l as t op l as t i c  analysis  when  the  peak  stress 
at  the  bole  exceeds  the  material  tensile  yield  strength. 


■-  --  ^  - i — -v. 


A  two-dimensional  finite  element  analysis  has  been  used  to  compute  the 
stress  distribution  in  the  uncracked  lug  for  lugs  of  various  Ro/Ri  ratios. 
The  major  useful  results  of  this  analysis  are  the  stress  distribution  along 
the  potential  crack  path  (Table  3-1)  and  the  stress  concentration  factor 
(Figure  3-5). 

The  compounding  method  combines  known  solutions  to  obtain  an  engineer¬ 
ing  approximation  for  the  stress  intensity  factor,  Equation  (10).  Stress 
intensity  factor  values  calculated  by  this  method  have  been  listed  in  Table 
3-2. 


The  two-dimensional  cracked  finite  element  method  properly  models  the 
crack  tip  stress  singularity  and  the  distribution  of  pin  bearing  pressure, 
which  changes  drastically  with  crack  length  (Figures  3-14  through  3-16). 

The  stress  intensity  factor  results  from  this  method  have  been  shown  in 
Figure  3-18  and  listed  in  Table  3-3. 

The  weighting  function  method  calculates  the  stress  intensity  factor 
as  the  integral  of  the  product  of  the  stress  in  the  uncracked  lug  times  the 
Green's  function  for  the  lug;  Equation  (14).  The  Green's  functions  for 
straight  lugs,  developed  using  two-dimensional  cracked  finite  element  analyses 
with  point  loads  applied  on  the  crack  surface,  have  been  listed  in  Tables  3-4 
through  3-6.  However,  these  "original"  Green’s  functions,  when  used  with 
the  stress  distribution  in  the  uncracked  lug  from  Table  3-1,  obtain  results 
at  variance  with  those  of  the  two-dimensional  cracked  finite  element  analysis. 
The  discrepancy  arises  because  the  Green's  function  method  fails  to  account 
lor  the  change  in  the  distribution  of  pin  bearing  pressure.  To  correct  the 
discrepancy,  the  original  Green's  functions  have  been  modified  such  that  the 
cracked  finite  element  results  are  exaetly  duplicated  when  these  'Modified" 
Green's  functions  (Tables  3-7  through  3-9)  are  used. 

Figure  3-29  compares  the  stress  intensity  factors  for  through-the-thiek- 
ness  cracks  compute  d  by  the  various  methods.  Reasonable  agreement  among  all 
methods  is  obtained  at  R0/Rj  ratios  of  2.25  and  3.0,  but  not  at  R^/R^  *  1.50. 
Thus,  any  method  could  be  used  lor  larger  R^/^  ratios,  but  the  sore  rigorous 
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two-dimensional  cracked  finite  element  method  (or  equivalently,  the  Modified 
Green's  function  method)  is  preferred  for  R^/R.  2.0. 

The  weighting  function  method  can  be  applied  to  account  for  residual 
stresses  caused  by  a  shrink-fit  bushing.  Assuming  the  lug  and  bushing  remain 
in  int imate  contact  during  loading,  the  stress  intensity  factor  is  calculated 
by  Equation  (24)  from  the  sum  of  the  residual  stress  caused  by  bushing  in¬ 
stallation  plus  the  distribution  of  stress  caused  by  the  applied  load.  The 
residual  stress  is  estimated  from  the  closed-form  solution  for  two  concentric 
cylinders.  Equation  (30).  The  applied  stress  is  obtained  from  two-dimensional 
finite  element  analysis  of  an  uncrackcd  lug  in  intimate  contact  with  a  neat- 
fit  bushing.  If  the  bushing  and  lug  are  of  the  same  material,  then  the  applied 
stress  distribution  can  be  obtained  from  Table  3-1.  Sample  applied  stress  dis¬ 
tributions  for  various  bushing/iug  modulus  ratios  are  shown  in  Figures  3-33 
and  3-34.  It  is  to  be  noted  here  that  this  method  can  give  unconservative 
results  when  separation  occurs  between  the  bushing  and  lug.  In  such  instances, 
an  improved  methodology  has  been  developed  and  is  presented  in  Volume  III  of 
t  h i s  repo r t  . 

"io  solution  methods  for  through-the-thickness  cracks  can  be  modified 
to  analyze  a  cornet  crack,  utilizing  a  corner  crack  correction  factor,  along 
with  a  method  to  account  for  the  transitional  behavior  as  the  corner  crack 
becomes  a  through-the-thickness  crack.  Two  alternative  correction  factor 
approaches  have  been  suggested,  a  one-parameter  and  a  two-parameter  method. 

In  the  one-parameter  method  the  tlaw  shape  is  assumed  to  be  constant 
(e.g.,  a/c  -  1.33)  and  the  stress  intensity  factor  at  the  lug  surface  point 
(Point  C)  is  calculated  using  Equations  (31)  and  (32).  These  equations  apply 
to  flu  crack  throughout  its  growth  from  a  corner  crack,  through  transition,  to 
a  through-the-thickness  crack. 

In  the  t wo -pa rame t e r  method.  Equal  ions  (  (!)  through  (36)  (in  conjunct  ion 
with  Figures  1/3,  l 76  and  1 7 8  ol  Ret c rente  [16])  are  used  to  compute  stress 
Intensity  factors  at  the  lug  surface  and  hole  wall  (Points  G  and  A).  During 
the  transition  to  a  t brought  hr -l hickness  crai k,  Equations  (37)  through  (40) 


arc  used  to  compute  stress  intensity  factors  at  the  front  and  back  surface. 


A  three-dimensional  cracked  finite  element  method,  although  too  ex¬ 
pensive  for  general  application,  has  been  used  to  check  the  accuracy  of  the 
one  and  two-parameter  corner  crack  solutions.  Comparisons  of  stress  intensity 

factors  are  shown  in  Figures  3-46  throuch  3-50  and  comparisons  of  corner  crack 
correction  factors  are  shown  in  Figure  3-51. 

An  elastoplast ic  analysis  has  been  described  for  use  when  the  peak 

stress  in  the  uncracked  lug  exceeds  the  material  tensile  yield  strength. 

Using  the  stress-strain  curves  of  Figure  3-52,  an  iterative  finite  element 

analysis  with  incremental  loading  and  unloading  is  used  to  calculate  the 

stress  distributions  in  the  uncracked  lug  for  the  maximum  and  minimum  loads 

of  the  fatigue  cycle.  These  stress  distributions  are  used  with  the  modified 

Creen's  function  to  estimate  K  and  K  ,  for  a  lug  with  a  throueh-the- 

max  min  ° 

thickness  crack.  Sample  results  for  AK  have  been  given  in  Figures  3-59  and 

3- 60  and  Table  3-13.  This  nonlinear  method  is  inexact  and  only  an  approxima¬ 
tion,  because  strictly  speaking  the  validity  of  the  Green's  function  method 
requires  linearity  between  load  and  stress. 

2,  SUMMARY  OF  STRESS  INTENSITY  FACTORS  FOR  TAPERED  LUGS  (SECTION  TV) 

Section  IV  covers  the  calculation  of  stress  intensity  factors  for  ta¬ 
pered  attachment  lugs  subjected  to  either  axial  or  off-axis  loading.  Only 
unbushed  lugs  with  a  45-degree  included  taper  angle  are  analyzed. 

A  two-dimensional  finite  element  analysis  has  been  used  to  calculate 
r he  stress  distribution  in  an  uncracked  lug.  For  axial  loading,  the  peak 
tangential  stress  is  located  at  approximately  90  degrees  to  the  lug  axis 
(Figure  4-5),  and  the  stresses  along  the  90-degree  line  are  listed  in  Table 

4- 1  tor  five  values  ot  R  /R  .  Kquat'on  (45)  is  an  equation  for  stress  con¬ 
centration  factor  tor  axial  loading  of  a  tapered  lug  with  a  taper  angle  be¬ 
tween  0  and  45  degrees.  For  off-axis  as  well  as  axial  loaOing,  the  angular 
locations  of  peak  stresses  are  given  in  Figure  4-5. 

Stress  intensity  factors  for  tapered  lugs  have  been  calculated  by  two- 
tit  mens  j  ana  1  cracked  finite  element  analysis.  Stress  intensity  factor  results 


162 


are  listed  in  Tnble  4-2  for  axial  tension  loading  (and  plotted  in  Figure  4-20 
for  axial  compressive  loading)  of  tapered  lugs  having  R^/R,  ratios  from  1 .  r> 
to  3.0.  Off-axis  loading  solutions  for  stress  intensity  factors  are  plotted, 
in  Figures  4-22  and  4-23  for  a  tapered  lug  with  R^/R,  -  2.25,  loaded  in  the 
-45  and  -90  degree  directions,  with  a  crack  at  either  of  the  two  most  criti¬ 
cal  locations. 

For  corner  cracks  in  tapered  lugs,  the  corner  crack  correction  factors 
presented  in  Section  III  can  be  used  in  conjunction  with  the  appropriate 
through-the-th ickness  stress  intensity  factors.  For  problems  involving  bush¬ 
ings  or  stresses  above  the  material  yield  strength,  the  weighting  function 
would  be  needed  for  tapered  lugs  for  each  new  loading  direction  and  crack 
orientation.  As  an  approximation,  the  Green's  function  for  the  axially- 
lo. uleil  straight  lug  with  the  same  R  /R.  ratio  may  be  used  although  the  ac¬ 
curacy  of  this  approximation  is  questionable,  particularly  for  off-axis  load¬ 
ing. 

3,  SUMMA R Y  OF  FATIGUE  CRACK  GROWTH  ANALYSIS  METHODS  (SF.CTTON  V) 

Crack  growth  analysis  methodology  has  been  briefly  summarized  in  Sec¬ 
tion  V.  Alternative  constant  amplitude  fatigue  crack  growth  rate  equations 
described  include  those  of  Paris  [49],  Forman  [50],  and  Walker  [51].  Alter¬ 
native  spectrum  load  interaction  models  summarized  are  the  Wheeler  [52], 
Willenborg  [53],  Generalized  Willenborg  [54]  and  lisu  [55]  models.  The  selec- 
t  ion  of  a  rate  equation  and  retardation  model  is  left  to  the  discretion  of 
the  analyst  who  is  conducting  the  era- k  growth  analysis. 
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APPENDIX  A 


MODE  I  STRESS  INTENSITY  FORMULAS 
RELATED  TO  CRACKED  LUCS 


LINE  CRACKS  IN  OOUBIY  INFINITE  PLATES 


CASE  11.  Ur’ifor'n  T*njil*  Strets 


CASE  12.  Centrai  Splitting  Porcil  on  Cftd»  Surtactt 


whan  t  •  thicknau 


Numncal  Esampla:  »0*  3,  P  ■  t,t  ■  0.6'  *  0.54289 


CASE  13.  Arp.t'irv  Splitting  Fprcil  pn  Cf»{k  Surf*C« 


NOTE.  *  ii  ntgativt  if  fore*  n  applied 
closer  to  track  tip  (1). 


(1) 


121 


Numerical  Esimpit:  a  •  3.  *  •  1,P  •  l.t  •  0.8:  K,j  •  0.76778,  K^j  •  0.38338 
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LINE  CRACKS  IN  INFINITE  LENGTH  STRIPS 


CASE  21.  Symmetric  Crack,  Uniform  Tension 


K21  ?  K11  *21  1 2b  ° 

$21(V  J/sec(V 


Numerical  E nample. 

a  «  3,  b  =  8.S  *  !.  K,.  *  3.3668 

0  ii 


CASE  22.  Symmetric  Crack.  Central  Spl  ttmg  Forces  on  Crack  Surfaces 


k22  *  Ki2  *21 

U22(V  Im^T 


Numerical  Example 

sQ  -  3,  o  *  3.  P  »  l.t  *  0.6  K22  *  0.63125 


CASE  35.  Single  Crack.  Distributer!  Splitting  Forces  at  Hole 


(1) 


K35  =  K13  *31  (P1*  u35  (P1* 


u35  (P,)  *  exp  0.15  IP,  -  1) 


x  =  aj/2,aQ  5  r  *  x 


limitation:  Total  load  P  results  from  uniform  radial  pressure  from 
tr/8  to  7ir/S  and  ir/8  to  -7ff/8.  (K  drops  drastically  if  pressure  is 
applied  between  -rr/8  and  rr/8). 


CASE  26.  Unsymmetric  Crack,  Uniform  Tension  (bn<.  bf) 


t  t  t  tst  t 


|-*-ao— *■- 

»-aO-*-J 

1  (f) 

(n) 

—  ^ 

1  -  U_  .  ^ 

1 

114  4S4  4 


At  (f): 


c(,) 

26 


11 


Attn):  K 


(n) 

26 


K  <t(n) 
11  26 


^Xf,  Xn^ 


ysec(Xnf)  -1 

- 

1  +0.21  sin  /  8  arctan 

r  "  _  x  o.9-| 

/  X  n  X  f  \ 

L  \ 

K+v  J 

ij 

CASE  27.  Unsymmetric.  Crack,  Central  Splitting  Forces  on  Crack  Surfaces 


K!7  '  K t3  2G  <X| '  X  1,1  ”z2‘Xnl 


K®  ■  k'2,1  4  "gI\,.Xn>“I2  'X»|t 


x  “  ^  ^n  Combining  of  Cases  12,  13,  21 

22  and  26. 

j  1  +  i>  :  P.  Stress 
^  P.  Strain 
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CASE  44.  Single  Crack.  Distributed  Splitting  Forest  at  Hole 


(1) 


K44  “  K27<t'31(P1l4,21!Vu3StP1) 


Limitation:  Same  load 

Combining  ot  Casas  13.  21. 

distribution  at 

27.  31  and  35 

Casa  3b. 

Dependent  Variable!: 


a  3  a,/2.  bB  3  b  -  a  and  aQ  3  r  ♦  a^/2 


CORNER  CRACKS,  UNIFORM  TENSION 


CASE  71,  Quarter-Circular  Corner  Crack  Which  Becomes  Through  Thickness 


t  t  s  f 


<C) 

71 


*<C) 

71 


0-2386 
1  ♦  2  la/()2 


NOTE  F  or  cracks  beginning  as  quarter  circular  corner  cracks  in 
configurations  32  through  70. multiply  the 
through-crack  solution  by^'^l 


Limitations: 

No  plate  bending 
a  =  c 


Numerical  Example  (0) 

a  =  c  2.  r  -  1,  t  =  3,  S  =  1:  Kn  =  2.18834 
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